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1 Introduction to the Problem 


1.1 Purpose of the Document 

This document is a product of a research project initiated in February 1999 by the X-38 Flight 
Controls Branch at the NASA Johnson Space Center (JSC). Funded by NASA Grant NAG9- 1 085, 
the effort was associated with the Flight Mechanics Laboratory (FML) of the Texas Engineering 
Experiment Station - the research arm of the Dwight Look College of Engineering at Texas A&M 
University (TAMU). One of the tasks of the unsolicited proposal that led to this grant was to provide a 
set of design guidelines that could be used in future by JSC. The subject of these guidelines was to be a 
flight control design for vehicles operating across a broad flight regime and with highly nonlinear physical 
descriptions of motion. The guidelines specifically were to address the need for reentry vehicles that 
could operate, as the X-38 does, through reentry from space to controlled touchdown on the Earth s 
surface. The latter part of controlled descent was to be achieved by parachute or paraglider - or by an 
automatic or a human- control led landing similar to that of the space shuttle Orbiter. 

Since these guidelines address the specific needs of human-carrying (but not necessarily piloted) 
reentry vehicles, they deal with highly nonlinear equations of motion, and their generated control systems 
must be robust across a very wide range of physics. Thus, this first- generation document deals almost 
exclusively with some form of dynamic inversion (Dl), a technique that has been widely studied and 
applied within the past 25 to 30 years. Comprehensive and rigorous proofs now exist for transforming 
a nonlinear system into an equivalent linear system. (Called either feedback linearization or Dl, it is 
based on the early papers of Krener and Brockett 1 2 .) At about the same time, theoretical advances 
essentially completed the background for ensuring the feedback control laws that make prescribed 
outputs independent of important classes of inputs; namely, disturbances and decoupled control 
effectors. These two vital aspects of control theory - noninteracting control laws and the transformation 
of nonlinear systems into equivalent linear systems — are embodied in what is often called Dl. Falb and 
Wolovich’ considered noninteractions as a facet of linear systems theory. Singh, Hugh. Freund, and 
Porter 4 ’ 5 h extended these notions into nonlinear systems. Isidori and his colleagues ' contributed 
significantly to Dl theory by using mathematical notions from differential geometry. Balas and his 
colleagues applied these ideas to a variety of aerospace flight control system designs - including the F- 
18 high angle-of- attack research vehicle (HARV) 9 as well as to the X-38"’ itself. They also provided 
powerful, commercially available software tools' 1 2 that are widely used by control design practitioners. 
Though there is no doubt that the mathematical tools and underlying theory are available to industry and 
government agencies, there are open issues as to the practicality of using Dl as the only (or even the 
primary) design approach for reentry vehicles. Our purpose, therefore, is to provide a set of guidelines 
that can be used to determine the practical usefulness of the technique. 

This document will answer the following questions related to four main topics: 


1 . If we use Dl as our primary design method, what tools are available to implement the design 
tasks? 

2. How easy is it to obtain and to learn to use these tools? Can an entry-level (an 
undergraduate) engineer be expected to be familiar enough with the tools to be productive 
without receiving specialized training and consulting help? 




3 . Is it easy to convey the value of using DI? How does a design group communicate the 

validation of systems modeled with this modem control technique? 

4. What form of robustness analysis is appropriate? Is more than one technique worth 
considering? 


Section 2 of this report addresses the first question by first summarizing the value of three tools 
used by TAMU FML engineers - MACH [Mutli- Application Control], MATLAB, and batch 
simulations. This section goes on to investigate and explore the available forms of robustness analysis 
(question 4) as the forms relate to practical uncertainties and disturbances. Section 2 concludes with 
first thoughts on how we would go about evaluating the various tools. 

Section 3 addresses how Dl is achieved from the perspective of new graduate students who has 
to teach themselves these techniques. It is hoped that later studies will expand and extend this validation 
process to show that less- sophisticated talent can also successfully complete workable designs. 

Section 4 illustrates the simulation component buildup surrounding Dl, and it applies Dl to the 
X-38 reentry vehicle model in three separate examples. The first tests a DI controller against a 
nonlinear MATLAB simulation to evaluate performance; the second and third present longitudinal and 
lateral/directional DI controller designs, respectively. 

Section 5 describes two different controller analysis techniques and analyzes DI controllers 
using both methods. The controller analysis techniques addressed in this section include |i- analysis and 
linear quadratic performance index analysis. 

Section 6 provides a summary of the theoretical background needed to understand some of the 
DI design procedures and to complete at least elementary robustness analyses of the DI system. 

Finally, Section 7 is a fairly extensive list of references used to prepare this report. Although the 
bibliography is not comprehensive, it does include much of the classical work that has been done to this 
point. 
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2 Synthesis Procedure 

Synthesis is the process by which the components or elements of a system are brought together 
by a designer to accomplish the tasks under consideration. The trick is to be sure that the individual 
parts are integrated in such a way that the sum of the parts produces an outcome greater than the 
individual contributions of the parts. This “synergy” is a result of an integrated design. Integration 
begins w'ith the process used, and depends strongly on the tools available. 

in this first iteration of our design guidelines, we will consider three sets of software tools; i.e., 
MACH, MATLAB, and batch simulations. MATLAB is a widely used commercial software package 
for control system design that has both a command line and a graphical user interface. It is also 
relatively easy to use, and many colleges and universities teach undergraduate courses that integrate 
MATLAB-based problems into their pedagogy. Moreover, MATLAB has a number of specialized 
toolkits that directly address matrix algebra and modem control system design, including D1 and 
techniques often used to analyze the robustness of such designs. This set of tools is quite extensively 
documented; indeed, MATLAB has steadily evolved and been improved over several years of 
commercial usage. 

MACH is a set of proprietary software tools developed and used (but not sold commercially) 
by Honeywell that directly address some of issues common to DI. One of the key questions we w ant to 
answer in this report is: Is it feasible for a relative beginner to build up DI models without using tools 
such as MACH? Or, is MACH indispensable to the efficient generation of DI modules? 

Finally batch simulation, which can be done (at least partially) within MATLAB’s Simulink 
module, is a software tool that requires some attention. It is doubtful that a control system designer 
today would attempt to produce a flight -worthy system without first generating at least a mathematical 
model of the specific system under consideration. We used the shuttle engineering simulator (SES) as 
the basis for our batch simulation buildup. As is almost always the case, keeping the simulation current 
as the vehicle (in this case, the X-38) design evolves is a recurring headache. As new data become 
available, the simulation has to be updated. Lesson Learned 1: Set up a procedure early in the 
process for updating and formatting aerodynamic (and other) databases. A corollary to this is that 
time and resources must be devoted to maintaining these databases or all facets of the program will 
suffer. Flight control design cannot proceed efficiently without this effort. 


2.1 Tools 

2.1.1 MATLAB 

As mentioned earlier, MATLAB is one of the most widely used commercial software packages 
available for control system design. So, it and its companion product, Simulink, are used extensively in 
this study. The DI controllers are developed in MATLAB, and simulations are run within the Simulink 
environment. This simulation development process and its results are at the core of this study. Detailed 
linear robustness analyses of the example controller for the X-38 are also made easier by MATLAB. 
Two different MATLAB toolboxes will facilitate p-analysis: the (1 ) Robust Control Toolbox and (2) u- 
Synthesis and Analysis Toolbox. Obviously, the latter is intended for p-analysis since it is dedicated to 
that process. Hence, a number of useful functions are readily available and packaged with a detailed 
instruction manual. A detailed discussion of the underlying theory of p- analysis is given in Section 5.1. 

2.1.2 Multi-Applica tion Control 

One of the spacecraft controllers using the DI approach is based on MACH. MACH is 
a proprietary software package developed by Honeywell that was previously applied to several flight 
control designs such as the F- 1 8 HARV and the X-29 aircraft. Its basic structure consists of an inner- 
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loop DI controller wrapped around an outer- loop classical proportional integral (PI) controller. Figure 
2.1 below shows the similarity between the MACH system structure and the controller used in the first 
example (see Section 4.6). Although slight differences exist, these are associated primarily with different 
definitions for the control variables (CVs). 



Figure 2. 1 Typical MA CH system structure 


Since MACH is proprietary to Honeywell, implementation details cannot be presented in 
a document of unlimited distribution such as these guidelines. Details are sketchy in any event; and other 
than a comprehensive outline of the MACH structure, the code is not used broadly as a design tool in 
this version of our guidelines. Later implementations may include more on the MACH software if it 
becomes obvious that the code is useful to the overall design process. For purposes of this document, 
we will therefore focus on demonstrating that Dl can be successively implemented with other tools and 
procedures with only a moderately intense learning effort on the part of the analyst. A brief comparison 
of a MACH controller and a Simulink example are presented in Section 4.4. This comparison primarily 
highlights the differences. 

2.1.3 Batch Simulation 

The batch mode of SES X-38-V201 version 1 .3 is used in this study. This version contains the 
shuttle-derived classical controller designed at JSC by John Ruppert. Although different versions of 
MACH have already been implemented for the X- 3 8- VI 32; this is the only version already 
implemented that uses the V201 database prior to release of the MACH controller in late 1999. This 
SES version was thus our only available choice since the scope of our study was to examine the 
characteristics of a DI controller throughout the entire X-38 flight envelope (i.e., from hypersonic 
through subsonic flight regimes of anticipated trajectories). 

By using batch implementation, a nominal case is executed to obtain essential vehicle properties 
(aerodynamic coefficients, moment/product of inertia, etc.). The data thus obtained are then 
incorporated into the Simulink- based DI controllers. A few attempts are also included in which we 
began examining uncertainty in the mass properties of the X-38 by varying these parameters slightly 
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(Section 4.6). However, time constraints as well as the complexity of the SES batch simulation limited 
the number we perfonned of these runs. 

2.2 Specifications 

The controller design procedure is normally iterative and centers around designing a controller 
that satisfies a set of design specifications. These specifications can be provided in the time domain, the 
frequency domain, or both. It is important to note that the type of input (i.e., impulse, step, ramp, 
sinusoid) must be specified. Examples of these are introduced and discussed below. 

2.2.1 Time Domain 

Time domain inputs, such as a step input, can be used to evaluate system characteristics such as 
damping, natural frequency, overshoot, etc. Initial condition or impulse excitations are particularly useful 
in evaluating the damping of rate variables. Such time domain controller responses can be evaluated 
through simulation in the Simulink environment, for example. 

2.2.2 Frequency Domain 

The response of a linear system to a sinusoidal input is referred to as the system s frequency 
response. Frequency domain specifications are concerned with the response of a system to frequency 
varying inputs, most often of the sinusoidal type. Typical specifications are gain margin, phase margin, 
and bandwidth. Gain margin is the amount by which system gain can increase before the system 
becomes neutrally stable. Phase margin is the amount by which phase lag can increase before the 
system becomes neutrally stable. Bandwidth - defined as the maximum frequency at which system 
output will satisfactorily track a sinusoid input - is basically a frequency domain measure of response 
speed. It is therefore akin to the time-domain specification of rise time. Frequency domain 
specifications are important to multiple- input, multiple-output (MIMO) robust controller design since 
most available methods are based in the frequency domain and thus use some or all of the frequency 
domain specifications. 

2.3 Uncertainty Modeling 

Actual controllers are expected to perfomi well for an entire class of transfer functions 
representing the range of plant dynamics and operating environment. Since it is impossible to 
analytically or empirically model with 100% accuracy a dynamic system and the effects of its operating 
environment, uncertainty modeling plays an important role in controller design and analysis. But even 
when applying optimal control techniques, the resulting controller designs will not be truly “optimal’' 
because 

• operating environments can introduce undesirable/unknown performance. 

• the system is inherently nonlinear (i EXAMPLE : Coulomb friction, hysteresis, backlash, and 
deadbands). 

• physical components are subject to wear and failure. 

• there are limitations to implementation ( EXAMPLE : computational delays). 

Uncertainties are broadly classified in two categories - structured and unstructured - both of w hich are 
usually present in any given physical system. The key to successful uncertainty modeling (and, thus, to 
robust controller design) is to recognize to which category a particular type of uncertainty belongs and 
then to determine the characteristics of that uncertainty. Further, the majority of robust control 
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techniques require uncertainty be modeled entirely in the frequency domain. These topics are outlined in 
the following sections. 


2.3.1 Structured 

Structured uncertainties can be modeled and have relatively well-known bounds and ranges. 
Parametric uncertainties arise from and include 

• control effectiveness, 

• aerodynamics, 

• mass, and 

• inertia. 


Prediction Accuracy of Aircraft Stability Derivatives Prediction Accuracyof Spacecraft Inertias 
Lift Curve Slope C L ± 5% 1 un ± 10% 

Pitch Damping C m ± 20% I„ p + 10% 

Yawing Moment Due To Roll Rate C n ± 90% 

Parametric uncertainties are important; but since by definition they can be understood and modeled, 
they can usually be analyzed and handled. Nonparametric uncertainties are potentially more dangerous 
because they are not as well understood and are difficult or impossible to model accurately. However, 
those that can be represented as some type of stochastic process can be easily incorporated into the 
design model. The key issue is to determine the relative magnitude of nonparametric uncertainties. In 
general, small nonparametric uncertainties cause small errors while large nonparametric uncertainties 
cause large errors. It is also important to determine how quickly nonparametric uncertainties vary. 

2.3.2 Unstructured 

Unstructured uncertainties are those for which generally little to no knowledge is possessed. 

They are usually not modeling- related nor can they be modeled at all. Nonparametric unstructured 
uncertainties include 

• high- frequency unmodeled dynamics, 

• actuator dynamics, 

• structural vibrations, 

• measurement noise, 

• round-off error and truncation, and 

• sampling delay. 

Since retention of full, nonlinear dynamics severely restricts the number of synthesis techniques 
presently available, linearization of actually nonlinear dynamics is often required. So, approximations are 
inherent and introduce uncertainty. Actuators also fall into this category for the same reason. 

Uncertainty due to structural vibrations and measurement noise can be represented with a certain degree 
of accuracy when experimental data is available. In the absence of experimental data or when a simpler 
representation is wanted, measurement noise is often approximated as a sine wave. Round-off error 
and truncation are extremely difficult means of representing uncertainties. No widely accepted standard 
method exists for them. 
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2.3.3 Frequency Domain 

Classical Control addresses the issue of uncertainty by assuming that all types of uncertainties in 
the system cause only gain changes, or phase changes, to occur. Robust Modem Control takes a 
frequency domain approach using transfer functions in the S- domain such that certain types of modeling 
errors are assumed to have certain frequency effects. Since parametric modeling errors are structured 
uncertainties with known bounds, they are assumed to cause low-frequency effects. Consequently, 
neglected and possibly higher- order dynamics are assumed to cause high-frequency effects. 

Unstructured uncertainties, which are not well understood, represent systems in the frequency domain 
whose frequencies simply are assumed to lie between some upper and lower bound. Additive 
uncertainty is used to model errors in neglected high- frequency dynamics; this represents the absolute 
error in the model. Multiplicative uncertainty, which is used to model errors in actuators or sensor 
dynamics, represents the relative error in a model. This latter type of uncertainty is most useful in 
robustness analysis and design. 

2.4 Disturbances 

Disturbance rejection properties to exogenous disturbances - e.g., gusts, turbulence, wind shear 
- are particularly critical in flight control system design. By definition, an exogenous input is one that a 
controller cannot manipulate. These unstructured uncertainties are stochastic processes and, as such, 
are best represented as stochastic models in terms of mean and variance. The standard gust and 
turbulence models, due to Von Karman and Dryden, are empirically based and directly applicable to 
both controller design and controller analysis. 

2.5 Dynamic Inversion Synthesis 

D1 synthesis is a controller synthesis technique by which existing deficient, or undesirable, 
dynamics are canceled out and replaced by desirable dynamics. Cancellation and replacement are 
achieved through careful algebraic selection of the feedback function. For this reason, this methodology' 
is also called feedback linearization. It applies to both single-input, single-output (S1SO) and MIMO 
systems, provided the control effectiveness function (in the SISO case) or the control influence matrix 
(in the MIMO case) is invertible. The method works for both full-state feedback (input-state feedback 
linearization) and output feedback (input-output feedback linearization). A fundamental assumption in 
this methodology is that plant dynamics are perfectly modeled and can be canceled exactly. In practice 
this assumption is not realistic, so the new dynamics require some form of robust controller (see Section 
2.6. 1 ) to suppress undesired behavior due to plant uncertainties. Examples of D1 synthesis are shown in 
Chapter 3. 

2.6 Robustness 

Compensators are designed to satisfy specified requirements for steady- state error, transient 
response, stability margins, or closed- loop pole locations. Meeting all objectives is usually difficult 
because of the various tradeoffs that have to be made and because of the limitations of design 
techniques. For example, although classical root locus design places a pair of complex conjugate poles 
to meet transient response specifications, the designer has little control over the location of all other 
poles and zeros. The particular property that a control system must have to operate properly in realistic 
situations is called robustness. A control system that possesses both good disturbance rejection and 
low sensitivity is said to be robust. Disturbance rejection is the ability to maintain good regulation 
(tracking) in the presence of disturbance signals. Low sensitivity is the ability to maintain good 
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regulation (tracking ) in the presence of changes in plant parameters. Mathematically, this means that a 
controller must operate satisfactorily for not just one plant but for a family or a set of plants. 

Robustness is divided into two distinct yet related categories: stability’ robustness and 
performance robustness. Stability robustness is the ability to guarantee closed- loop stability in spite of 
parameter variations and high-frequency unmodeled dynamics. It is important to note that relative 
stability, not absolute stability, is of interest in this context. Performance robustness is the ability to 
guarantee acceptable performance (settling time, overshoot, etc.) even although the system may be 
subject to disturbances. The Classical Control method quantifies robustness through gain margin and 
phase margin. Modem Control techniques use the structured singular value analysis of Section 6.2.3 to 
quantify robustness. In the MIMO case, both the maximum and the minimum singular values are 
measures of the amplification and attenuation, respectively, of the transfer function matrices that 
represent the family or set of plants of a system. Section 5.1 presents this robustness technique and 
demonstrates how to perform the analysis and interpret the results. 

2.6.1 p-Synthesis and H„ 

Structured singular value synthesis, or )U- synthesis, is a multivariable design method that can be 
used to directly optimize robust performance. It involves both p- analysis and FL synthesis. 

Performance specifications are weighted transfer functions describing the magnitude and frequency 
content of control inputs, exogenous inputs, sensor noise, tracking errors, actuator activity, and flying 
qualities. A family of models (consisting of a nominal model plus structured perturbation models) is used 
with magnitude bounds and frequency content specified using weighted transfer functions. All of this is 
wrapped into a single standard interconnection structure that is then operated upon by the algorithm. 

The PL control controller design methodology is a frequency domain optimization for robust 
control systems. FL is defined as the space of proper and stable transfer functions - i.e., transfer 
functions with a number of zeros less than or equal to the number of poles. The objective is to minimize 
the PL norm. Physically, this corresponds to minimizing the peak value in the Bode magnitude plot of 
the transfer function in the SISO case or the singular value plot in the MIMO case. There are certain 
advantages in minimizing the infinity- norm. These are 

• The infinity- norm is the energy gain of the system. By comparison, the Linear Quadratic 
Gaussian (LQG) technique minimizes the 2-norm, which is not a gain. 

• The infinity-norm minimizes the worst-case root mean square (RMS) value of the regulated 
variables when the disturbances have unknown spectra. The 2-norm minimizes the RMS 
values of the regulated variables when the disturbances are unit- intensity, white noise 
processes. 

• PL control results is guaranteed stability margins (and is therefore robust), whereas LQG 
has no guaranteed margins. 

As in the Linear Quadratic Regulator (LQR)/LQG methodology, PL is iterative. In the standard 
problem, the solution for the infinity-norm is iterated upon until it is less than a specified scalar value, 
gamma - known as the gamma iteration. In the optimal problem, the infinity- norm is progressively 
reduced until a solution does not exist. In the FL control problem, the weights are the only design 
parameters the user must specify. Constant weights are used for scaling inputs and outputs. Transfer 
function weights are used to shape the various measures of performance in the frequency domain; 
weights are also used to satisfy the rank conditions. Proper selection of weights depends a great deal 
on understanding both the modeling process and the physics of the problem. 

Necessary conditions for a solution are the ability to stabilize and detect the system; to perform 
various rank requirements on system matrices; and to ensure that the transfer function between 
exogenous system inputs and the outputs remains nonzero at high frequencies. This last condition, which 
is often violated, occurs because the transfer function is strictly proper; i.e., has more poles than zeros. 
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Solutions to H. and LQG problems are very similar. Both use a state estimator and feed back the 
estimated states, and both solve two Ricatti equations to compute controller and estimator gains. The 
difference in the solutions lies in the coefficients of the Ricatti equation and in an extra term in the 1 f 
solution. Examples of this methodology are presented in Chapter 5. 

2.7 Validation 

Validation - which consists of an attempt to match outputs between two different control and 
simulation software packages for the same control inputs, and for the same controller structure and gains 
— was performed on all examples in this document to ensure as much fidelity as reasonably possible. 

The degree of fidelity depends on the purpose of the example, the software tool used to synthesize and 
simulate the example, the operating system and language, and the platform on which the example was 
being run. 


2 . 7 . 1 MA TLAB versus MACH 

MATLAB and MAC H have similar structures that, in theory, should permit good validation. 
MATLAB was run on a personal computer (PC) and MACH was run on a UNIX-based workstation. 
The difficulty involved with this validation effort stemmed from a lack of understanding of the MACH 
code itself due to a lack of documentation. Although agreement between the two codes was generally 
good, it was inadequate for in-depth investigations and research. 

2 . 7. 2 MA TLA B versus Ba tch Simulation 

MATLAB was run on a PC, and the SES batch simulation was run on a Silicon Graphics 
Incorporated (SGI) UNIX workstation. Because of adequate documentation and open access to the 
SES source code, validation between these two codes proceeded rapidly and with excellent agreement. 
These two software codes form the basis for all of the controller design research presented in this 
document. 
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3 Applying Dynamic Inversion 


3.1 Introduction and Philosophical Approach 

This section shows how DI is applied to a relatively simple aircraft control problem. As will be 
explained in more detail in Section 4, since the concept of DI itself is quite simple, a controller can be 
designed in many different ways. For example, the controller might be either linear or nonlinear. Also, a 
DI controller is not limited to a first-order inversion. It can take on higher-order forms as well. This 
chapter describes one way of designing a Dl-based controller. The steps taken in completing this 
design are carefully delineated in the hope that a step-by-step outline will help others design Dl-based 
controllers. 

First, a brief outline of the DI process will be given to quickly review the concept, followed by a 
detailed description of how to design each controller component. Then, aircraft equations of motion are 
introduced, and the DI design process is applied to a particular reentry vehicle; i.e., the X-38. Finally, 
several forms of desired dynamics are presented for this DI application. 

3.2 Dynamic Inversion Concept (Linear Aircraft Controller) 

As we suggested previously, the basic concept of DI is quite simple. In general, aircraft 
dynamics are expressed by 


j&=F(x,u) 
y = H(x) 


(3.1) 


where .v is the state vector, it is the control vector, and y is the output vector. For conventional uses 
(where small perturbations form trim conditions), the function F is linear in a. Equation (3. 1 ) can be 
rewritten as 


&=f(x) + g(x)u (3.2) 

where/is a nonlinear state dynamic function and g is a nonlinear control distribution function. If we 
assume g(x) is invertible for all values of x, the control law is obtained by subtracting /(x) from both 
sides of Equation (3.2 ) before multiplying both sides by g '(x). 

u=g-\x)[&-f(x)] (3.3) 

The next step is to command the aircraft to specified states. Instead of specifying the desired 
states directly, we will specify the rate of the desired states, $. By swapping # in the previous 
equation to J^ s , we get the final form of a DI control law. 

u= g 1 (x)[^ es -f(x)] (3.4) 

Figure 3.1 shows a block diagram representation of the DI process. 
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Figure 3. 1 Dynamic inversion process 


Although the basic Dl process is simple, a few points need to be emphasized. First, although 
we assume g(x) is invertible for all values of ,v . this assumption is not always true. For example, g (.\ ) is 
not generally invertible if there are more states than controls. Furthermore, even ifg(.v) is invertible (i.e., 
g(.v) is small), the control inputs, //, become large; and this growth is a concern because of actuator 
saturation. Since the dynamics of the actuators, as well as sensor noise in the feedback loop, are 
neglected during this primitive controller development to illustrate the process, a perfect inversion is 
not possible. 

Dl is also essentially a special case of model- following. While it is similar to other model- 
following controllers, a Dl controller requires exact knowledge of model dynamics to achieve good 
performance. Robustness issues therefore play a significant role during the design process. (This issue 
is discussed in detail in Chapter 5.) To overcome these difficulties, a Dl controller is normally used as 
an inner- loop controller in combination with an outer- loop controller designed using other control design 
techniques. 

The closed- loop transfer function for a desired CV that is being inverted is found according to 
Figure 3.2. From this block diagram, we can observe that the desired dynamics operate on the error 
between the commanded CV and its feedback term. In this figure, the pure integrator on the right side 
is used to approximate the rest of the system dynamics, as shown on the right side of the block 
diagram 12 . The CV here corresponds to the state .v in the previous development as well as in Figure 
3.1. 


CV 


cmd j 




Desired 

Dynamics 


CV 


des 


r 
I 

I 
I 
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CV 


Actuators 
Dynamic Inversion 
Effector Allocation 
Airplane Dynamics 
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Figure 3.2 Block diagram to calculate closed-loop transfer function 


3. 2. 1 Simplified Longitudinal Controller for an Aircraft 

A simplified form of the linear longitudinal equation for an aircraft’s pitch axis considers only the 
pitching moment equation. 


<$= M„a+ M q q+ M Se 8e 


(3.5) 


1 1 









The correspondence of this equation to Equation (3.2) is clearly seen in the following mappings: 

x=> q 

u=>8 e 

f => M a a + M q q 

g=>M &e 

Since A/& is a constant for a linear time invariant system, the inverse of the control distribution function, 
g, is always obtained as a constant, 1/M&. 

Now, we need to invert this equation for the elevator deflection angle. Mapping is achieved by 
substituting the relationships shown in the correspondence mappings (above) into Equation (3.3), giving 
the following equation: 



(3.6) 


To obtain a control law, we specify the desired value of pitch acceleration, <$P S . Then, by substituting 
for (§■ in Equation (3.6) and by substituting a meas and q nieas for a and q, we get the following 
longitudinal dynamic inversion control law: 

gcmd = _ M^eas _ M ^ess j (3 .7) 

Wc ^ 

Figure 3.3 shows the block diagram representation of the longitudinal DI controller. 


qf des 



g cmd 


Figure 3.3 Longitudinal Dynamic Inversion Control block diagram 


Recall that aircraft dynamics are modeled as a simple first-order form (Equation (3.5)) 
to develop this simplified DI control equation. In this model, both nonlinearity and higher- order terms in 
the actual aircraft dynamics are neglected. Since this simple DI controller cannot completely cancel out 
the aircraft dynamics, controller performance is potentially degraded. 

Similarly, 8 e ^ 8 c e md due to actuator dynamics. This shortcoming, which is also neglected while 

simplifying the control law development, is most noticeable when the control surface position and rate 
exceed their limits - something that occurs often when the value of A/&. is too small (in this case, 

& c e md was unbounded). 
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Finally, a meas * a; q meas * q due to sensor processing. This factor is also neglected in the 
control law development, thereby potentially harming controller performance as well. 

3.2.2 Simplified Lateral Directional Controller for an Aircraft 

Lateral/directional DI control equations are developed in this section. Although the 
development procedure is similar to that of the longitudinal case, we need to simultaneously deal with 
two states (roll rate and yaw rate) controlled by two control surfaces (ailerons and rudders) instead of 
with one state (pitch rate) controlled by one control surface (elevator) as in the simplified longitudinal 
case. 

Simplified linear lateral aircraft equations can be written with respect to roll as well as yaw axes 
as 


/&= L p P + L r r + LpP + L 8a 5a + L^br 
&= N p p+ N r r + A/pp + N 6a ba + N^br 


If we write Equation (3.8) in a compact matrix form, we get 





o’ 






A 4 A 

l 

r 


4 a 4 r 

ba 

/& 


_N P N, A 

1 

A 


N 6a N dr 

br 


When we compare the matrix form of Equation (3.9) to Equation (3.2), each parameter is either a 
vector or a matrix but the form remains the same. 
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/& 

Pi 


x = 
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u - 


4 
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4 4 

N r N p 

4, " 

V 


(3.10) 


Notice here that the control distribution matrix, g, is a square matrix. Therefore, its inverse exists in 
general. 

As a next step similar to the longitudinal case, we will invert the roll rate and yaw rate dynamic 
equations to obtain aileron and rudder deflection angles. 
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(3.11) 
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Then, substituting the desired states /S^ es and for /# and $ along with the measured values of p, r, 
and (3 (p"‘ ea \ and (3 for p. r, and [3, we get the lateral DI control law. 
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(3.12) 


Figure 3.4 presents a block diagram representation of the lateral DI controller. 
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Figure 3.4 Lateral Dynamic Inversion Control block diagram 


3.3 Nonlinear Dynamic Inversion 

The previous examples illustrate DI control for a linear system. This approach can be readily 
extended to a system with nonlinear characteristics by starting with the following set of nonlinear 
equations typical for an aircraft. 


I 2 L+I X2 N UH+Qpq [/ 2 (/ y -l z )-l'j 2 ]qr 
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X* z xz 


Uz ~ r„ 


(3.13) 


Now, we will assume the longitudinal and lateral- directional moments - L, M, and jV- are linear with 
respect to aerodynamic derivatives; i.e., 
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(3.14) 


L = ^ + kb a +L Si S r +L pP +L r r 
M = M u a + M q q + M 6 b e 
N = A/p(5 + /V 6a 5 a + N K b r + N p p + N r r 


By substituting the above linear moment equations into Equation (3.13), we can obtain a relation in 
Equation (3.15) that combines linear and nonlinear tenns. 
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If the last term is ignored, tine result is identical to the linear set of DI equations previously obtained. 
Finally, inverting the above equation as well as performing proper substitutions of the commanded, 
desired, and measured values gives the resulting DI control law. 
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(3.16) 


3.4 Applying the Dynamic Inversion Controller to 
the X-38 - the Overall Structure 

The DI control laws developed in the previous sections are now integrated into an overall 
control structure. As the block diagram in Figure 3.5 shows, DI control is used as an inner loop 
accompanied by a and (t) feedback outer loops. Although any type of control technique can be used 
for the outer loop, simple feedback is used in this particular example to illustrate the characteristics of 
inner- loop DI control. 
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Figure 3.5 Overall Dynamic Inversion Control block diagram 


The overall DI controller requires commanded values of angle-of-attack, a‘"“\ and bank angle, 
<j) c "' rf , as inputs. Then, the measured values of a""’"' and b"""' are subtracted from the commanded 
values to produce a' "™ and in the outer loop. These error values are then fed into the Command 
Inverter block to be changed to rate commands, p‘""\ q‘"“\ and r 1 "" 1 . The Desired Dynamics block uses 
these rate commands and the rate measurements to create the desired acceleration terms - favored 
forms of commands for the DI controller. The next block is the DI block, which produces the control 
surface deflection angle commands 8/'" rf , 5/"" / , and 5/"" / . Finally, the control surface commands are 
fed into the Plant block, X-38 Model, via the Control Surface block. The Control Surface block 
includes control surface management logic, which blends the three command values, 8/"" / , and 
8<™" / , into two command values, 6 el "" 1 and 8,‘"' </ , that include the dynamics of the actuators as well as 
the position and rate limits of the actuators. Gust and sensor noises are added to the system as external 
disturbances as well. 

3.4. 1 Command Inverter 

In aircraft applications, sometimes it is better to command displacements in the angle-of-attack 
and bank angle rather than command the body axis rates p , q, and r. However, rate commands are 
needed as inputs to the Desired Dynamics block. The Command Inverter block (Figure 3.6) changes 
displacement commands into rate commands so that displacement commands are directly implemented 
in the DI controller. This section describes how displacement commands are transformed into rate 
commands. 
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Figure 3.6 Command Inverter block diagram 


3.4.2 Body Components and Euler Angles Relationship 

Roll, pitch, and yaw rates are obtained from Euler angular rates using the following 


transfonnation matrix: 


p I 1 0 -sinQ ^ 

q = 0 cosd) cosQsinq 

r 0 -sim . |> cos0cos<|) 


(3.17) 


Now, by substituting the commanded values p cml , if"'\ and for the corresponding p, q. and r and 
by replacing $, $, and with their corresponding commanded values, $ md , $ md , and 
\$£ md = 0 , the following relationship is obtained: 


' 1 0 
= 0 cos <j) 
0 -s/mb 


flQmd 


(3.18) 


The next step is to express commanded values of Euler rates in terms of the commanded values of the 
angle-of-attack and bank angles. 

3.4.3 Roll Angular Rate 

The commanded roll rate, <§£ md , is obtained from the commanded bank angle, o cmd . simply 
by differentiating with respect to time. 


p md _ foCmd 

dr 


(3.19) 


By substituting the above expression into the first row of Equation (3. 17), p"'" / is expressed as a function 
off""'. 

p cmd = d_.cmd (3.20) 

dt 
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3.4.4 Pitch Angular Rate 

Expressing pitch angular rate, , from angle -of- attack is slightly more complicated than the 
roll angular rate case. First, the Euler pitch angle can be expressed in terms of a (angle-of-attack), (3 
(sideslip angle), y (flight path angle), and 6 (bank angle) by 


0 = tan 


_f| ab + sin jyj a 2 - sin 2 y+ b 2 


sin 2 y 


(3.21) 


where- a = coso. cos (3 

b = sin <(> sin (3 + cos §sin a cos p 

The commanded value of the Euler pitch rate is calculated by differentiating the commanded value of 
Euler pitch angle by 


0 . 

dt 


(3.22) 


Substituting this expression for (f into the second row of Equation (3.17), t\'" d is expressed as a 
function of ff ml . 


qcmd 


= cosfy 


( cf q cmd 

[dt 


\ 

) 


(3.23) 


with 


0 ' 


cmd 


tan 



\ 


where: 


a cmcy =cosa cmd cosp 

b cmd = sim f md s/np + cos® cmd sina cmd cos p 


3.4.5 Yaw Angular Rate 

Instead of defining the corresponding Euler pitch and yaw rate commands to r 1 "" 1 , we simply set 
r "" 1 equal to zero. 


r cmd = 0 


(3.24) 
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3.5 Multiple Time Scale Method 

To bypass a singularity problem in the inversion of an ineffective control matrix, a multiple time 
scale method has been developed that has been found to be quite successful in solving the problem. 

This approach is especially useful when inverting slow- motion variables, such as angle -of- attack, a, in 
the longitudinal case and sideslip, (3, and bank angle, (J), in the lateral/directional case. These variables 
are deemed as “slow” dynamics because the control effectiveness on their dynamics is quite low. 
Variables making up the “fast” aircraft dynamics include pitch rate, q, in the longitudinal case, roll rate, 
p , and yaw rate, r, in the lateral/directional case. Since the control effectiveness on these body rates is 
high, these dynamics are considered “fast” dynamics. The multiple time scale method thus seeks to 
reformulate the original differential equation (Equation (3.1 )) into a set of two separate differential 
equations consisting of a set of slow dynamics, i&, and a set of fast dynamics, j&. 

&= f( x) + g(x)y (3.25) 

j&=h(x,y )+k(x,y)u (3.26) 


Applying this technique to the linear aircraft dynamics, i&= Ax + Bu , yields the following slow dynamic 
equations for the rate variables (Equation (3.27)) and fast dynamic equations fot the acceleration 
variables (Equation (3.28)): 
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(3.27) 


(3.28) 


where A andB represent the longitudinal state and control input matrix values for the linear state- space 
model, and A and B represent the lateral/directional state and control input matrix values. Also, the 
subscripts denote the row and column value, respectively. Note that in Equation (3.27), rate variables 
form the input for the slow dynamics while the actual control surface commands fonn inputs for the rate 
dynamics shown in Equation (3.28). Inverting each set of differential equations generates two DI 
control laws, one for the outer DI loop (Equation (3.29)) and one for the inner DI loop (Equation 
(3.30)). 
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(3.30) 
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A block diagram representation of this 2-time scale approach is shown in Figure 3.7. 


O^cmd Pemd Oa,cmd 



Figure 3. 7 Block diagram of the 2-time scale approach 


In the Fast Inversion block, fast desired dynamics are calculated and the control law in Equation 
(3.30) is implemented. Fast dynamics are a function of the CV commands, (p a „<i, Qamb and r cmil ) and 
their feedback terms (p, q, and r). Similarly, in the Slow Inversion block slow desired dynamics are 
calculated and the control law (Equation (3.29)) is implemented. Again, slow dynamics are a function 
of the CV commands (c i cmh (i and <| )„„,/) and their feedback terms (a, (3, and 6 ). In summary, the 
Slow Inversion block produces the commanded rate variables of Equation (3.29) that are fed to the 
desired dynamics in the Fast Inversion block. Using these fast desired dynamics, the fast inversion 
control law of Equation (3.30) produces the commanded control deflections that are sent to the control 
surface actuators, which then serve as input to the inherent dynamics. 

Several observations can be made from these two DI control laws. First, only the short- period 
aerodynamic terms (A 22 , A 23 , A 32 , and A 33 ) are present in this set of slow and fast dynamics. Further, 
these two equations combine, retaining all original lateral/directional state matrix terms. It is also 
important to observe that the control effectiveness of the elevon on angle-of-attack, B 2 i, is not present 
in the inversion matrix and has actually been eliminated altogether from these two sets of equations. This 
is the term that traditionally causes a singularity effect on inversion because the value is typically small in 
magnitude. Instead, the control effectiveness on the pitch rate dynamics, B 21 , has been retained for 
inversion in the fast DI control law. Similarly, control input matrix values affecting sideslip and bank 
angle dynamics have also been eliminated (B, N B ]2 , B 4 ,, and B 42 ). Therefore, only the control matrix 
terms for the rate dynamics have been kept (B 2 i, B 22 , B 3 |, and B 32 ). This is of benefit because the 
control surfaces are more effective on the rotational rate variables than they are on the rotational 
variables. Finally, it is important to emphasize the fact that this 2-time scale method requires that the 
designer specify two sets of desired dynamics: one set for the slow dynamics and one set for the fast 
dynamics. 
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3.6 Desired Dynamics 

The Desired Dynamics block, which was introduced during D1 control law development, is 
explained in detail in this section. 

D1 control requires acceleration terms. For example, as the following longitudinal Dl equation 
shows, a desired value of pitch angular acceleration, $, is required: 


cmd 

S e 



$ es -M a a meas 



(3.25) 


However, applications normally use either displacements or rates as command states to control the 
system. The Desired Dynamics block acts as a mapping function between the rate commands and the 
desired acceleration terms, which are the required form for the DI equations. The structure of the 
Desired Dynamics block is shown in the flow chart in Figure 3.8. 



Figure 3.8 Desired dynamics development for dynamic inversion 

(adapted from Ref. 12) 


Several forms of desired dynamics are presented in this document and are evaluated in terms of 
performance and robustness. The different forms of desired dynamics consist of 


21 







• Proportional dynamics 11 

• PI dynamics 12 

• Flying quality dynamics 1 4 

• Ride quality dynamics 

3. 6. 1 Proportional Case 

The simplest way of achieving desired dynamics implementation is the proportional, or first- 
order, case. In this case, the desired dynamics are expressed as 

CV*.-K„(CV_-CV). (3.31) 

The Ku term in Equation (3.3 1 ) sets the bandwidth of the response. The bandwidth must be selected 
to satisfy time-scale separation assumptions without exciting structural modes or becoming subject to 
the rate limiting of the control actuators. Figure 3.9 shows the block diagram representation of the 
Proportional Desired Dynamics block introduced in this section. 




Figure 3. 9 Proportional Desired Dynamics block diagram 



As shown above, the constant amplifies the error between the CV command and its feedback term. 
In Figure 3.9, CV is represented as the state, x. So, the closed-loop transfer function for the 
proportional form of desired dynamics, shown in Equation (3.32 ), desires to place a single pole at 


cv . K M 

CV cmd s + K w 


(3.32) 


3. 6. 2 Proportional In tegral Case 

The Desired Dynamics block is not limited to a first-order component. If the Desired Dynamics 
block does not create satisfactory handling qualities (for piloted aircraft) using a set of first-order 
equations, a higher-order system is used. A commonly used higher-order block is a PI. This form is 
particularly popular in DI literature that uses fighter aircraft examples 1215 . This type of Desired 
Dynamics block structure is also used in the linearized MACH controller designed by Honeywell for the 
X-38 vehicle and has been adopted for this study as well. The block diagram representation of a PI 
desired dynamics component is shown in Figure 3.10. It has the same form as that used in the 
Honeywell study 12 with a K R of 5 sec 1 selected. 
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The block diagram for the PI form corresponds to the desired dynamics in Equation (3.33 ), where x is 
the CV. 


CVdes=K, 


-CV cmd -CV 


'+'7 L ( CV --CV) 

4s 


These dynamics have a closed- loop transfer function of 


Ik, 


CV 


cv - s+ v 


(3.33) 


(3.34) 


which seeks to place a pole at s = — 0.5K,, for any real constant Kb. In essence, this form of dynamics 
compensates for both the CV and the CV rate. 

3.6.3 Flying Qualities Case 

Desired dynamics can also be specified in terms of flying quality levels. Mil- STD- 1 191 A u 
contains the flying quality specifications for different vehicle classes and mission types. Using this 
information, the proper time domain characteristics corresponding to a desired flying quality level 
(damping ratio, natural frequency, time constant) can be selected. These characteristics can be used to 
determine the proper values for the gains and poles locations in this form of desired dynamics. The 
flying qualities desired dynamics, shown in Figure 3.1 1, can be represented as 

CV„ s = K|q(S+a) (CV, u , -CV), (3.35) 

s' + bs + c 


where b = 2^ io w il dcs and c = 0 )^ dcs -K lq for the desired damping and natural frequency co nA . s . 
Both the gain, K t| , and zero location, a, are real constant values. 


23 




cmd 





Figure 3.11 Flying Qualities Desired Dynamics block diagram 


The closed- loop transfer function for the flying qualities dynamics is given in Equation (3.35). 
These dynamics desire to place three closed- loop poles and to add a single zero to the system. 


CV _ K r „(s+a) 

CV cmd s' , + bs 2 +(c+K t .)s+K 1 ,a 


(3.36) 


3.6.4 Ride Qualities Case 

The ride qualities form of desired dynamics that can also be used in D1 are given in Equation 
(3.37) and represented pictorially in Figure 3.12. 

CV*,=-Ei-(CV„.,-CV) (3.37) 

s+ b 


X 


cmd 




Figure 3. 12 Ride Qualities Desired Dynamics block diagram 


The closed- loop transfer function for this set of desired dynamics is given below by 

CV _ K rq 
CV cmd s 2 + bs + K rq 


(3.38) 


which desires two places to closed- loop poles at s, , = -0.5b ±0.5y b 2 - 4K ri| i for real constants b 
and K ri| . 

For highly augmented airplanes, the control anticipation parameter (CAP) replaces the 
longitudinal short period requirements, such as damping ratio and natural frequency' . The desired 
longitudinal dynamics are instead designed by selecting a desired damping ratio and CAP value. 
Equation (3.39) shows how to calculate CAP for a specific load factor n,, and natural frequency. 


CAP = — 


(3.39) 
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The flying quality levels for various CAP and i^r values are shown in Figure 3.13. 


Inherent 

Dynamics 



Figure 3.13 Control anticipation parameter requirements for highly augmented vehicle 


Once a CAP and a ^si> are selected to satisfy a desired flying quality level, the desired short- 
period natural frequency can be calculated from Equation (3.39). The gain and the pole location for the 
open- loop desired dynamics are then backed out from these two specifications, ^ S p and 0) n .sp. For 
example, the desired dynamics (see Figure 3.13) correspond to the following: 

• CAP = 0.802 

• Csp = 0.8 

Using the CAP equation yields a desired o) n , S p of 1.4s '. The corresponding gain and pole locations are 
then found according to Equations (3.40) and (3.41 ), respectively. 

K rq =0)i SP =1.96 (3.40) 

-b= -2^,^ = -2^ s ,=-2.24 (3.41) 

The desired dynamics for this example become 

c ‘ v *-=— L ^7< CV »- | - CV) - I- 1 - 42 * 

s + 2.24 
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3.7 Issues in Dynamic Inversion 


The procedure illustrating the main steps in DI controller design is shown below. Also listed 
next to each step are some of the issues that were encountered when first learning and using this design 
methodology. Some possible solutions, or options, to each of these issues are listed below the 
procedure. 


STEPS 


ISSUES 



Options: 

1 . If the inverse of the control input matrix does not exist, a multiple time scale method can be 
used — such as the 2-time scale approach presented in Section 3.5. A command inverter, 
such as those presented in Section 3.4, or a higher order of feedback linearization is also a 
possibility. 

2. There are no limitations on the form the desired dynamics may take. However, some of the 
common forms found in the literature include proportional, PI, and flying qualities. 

3. If redundant control effectors are available, a control allocation scheme can be designed in an 
effort to keep the required control deflections within the constraints of the actuator. 
Adjustment or replacement of the desired dynamics may also help reduce the control 
response. 

4. A robust outer loop is required because Dl alone does not guarantee robustness. The most 
popular robust outer- loop design methodology for Dl controllers is ji- synthesis. Although in 
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the literature LQG is not as well published in regard to DI controllers, LQG has been shown 
to be effective, as shown in Section 5.2. LQR is another possibility for robust outer- loop 
design as well. 
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4 Simulation 


In this chapter, the simulation component buildup is presented for the X-38 reentry vehicle, and 
the DI design process is then applied to the vehicle model. Simulation results are shown and discussed 
for three DI design examples to give the reader more detail on how the design process is carried out. 

4.1 Control Surfaces 

The Control Surfaces block (Figure 4.1 ) is composed of two sub-blocks, the Control Surface 
Management Logic sub-block and the Actuator Dynamics sub-block. A conventional aircraft is usually 
equipped with three control surfaces; namely, a rudder, an elevator, and an aileron. But, the X-38 has 
only two sets of control surfaces - rudders and elevons. Control Surface Management Logic must 
therefore blend three inputs into two inputs so that commands can be fed smoothly into the Actuator 
Dynamics sub-block. 



Figure 4. 1 Control Surfaces block diagram 


4. 1. 1 Definitions (from the X-38 Aerodynamic Design Data Book) 

The X-38 vehicle has two sets of control surfaces (see Figure 4.2): a pair of elevon control 
surfaces, located on the lower rear of the vehicle; and a pair of rudders, one at the top of each of the 
vertical fins. The positive direction of deflection for the elevon is down (as shown in Figure 4.2), and, 
looking from the rear, to the left for the rudders. Both elevons and rudders are dual- function control 
surfaces. Each surface is deflected independently to provide the required control authorities. 
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The elevon deflections are averaged to give the total elevon angle or elevator angle for pitch control. 

(8 e / _ +8 eR ) (4 i 


8 e = 


The average of the difference gives aileron angles for roll control. 

c. ( 8 e/_ — 8 ef?) 

s a = — 2 — 

Similarly, the rudders’ deflections are averaged to give total rudder for yaw control. 

( 8 r L +8 rf?) 

s f = — — 

The average of the difference gives speed brake angles for pitch and drag control. 

_(8 fjL -8 rR ) 

°sb = o 


(4.2) 


(4.3) 


(4.4) 


4.1.2 Control Surface Limits 

4. 1.2.1 Deflection position limit 

The X-38 control surface deflection limits are listed in Table 4.1. 
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Table 4.1 X-38 Control Surface Deflection Limits 


Bodyflap Lower Deflection Limit 

0 . 0 " 

Bodyflap Upper Deflection Limit 

45.0" 

Rudder Lower Deflection Limit 

-25.0” 

Rudder Upper Deflection Limit 

25.0" 


4. 1.2.2 Surface actuator rate limits 

Though the actual rate limits of the actuators used for the X-38 vehicle are set as a function of 
hinge moment, the constant values shown in Table 4.2 were used for this study. 


Table 4.2 X-38 Control Surface Rate Limits 


Body Flap Deflection Rate Limit 

50 deg/sec 

Rudder Deflection Rate Limit 

50 deg/sec 


4. 1.3 Control Actuator Modeling 

Control surface actuators are modeled with the following second-order lag for both rudders and 
body flaps: 


G, ictvator(s) ~ — 5 — p- (4.5) 

s z + 2^w n s + o)„ 

with co n = 26 rad/sec and C, = 0.707 in both sets of actuators. 

4.1.4 Control Surface Management 

Unlike conventional aircraft, the X-38 is equipped with only two pairs of control surfaces; 
namely, rudders and elevons. By commanding the deflections either symmetrically or asymmetrically, 
these two pairs of surfaces provide the same control effects that conventional rudders, elevators, 
ailerons, and speed brakes provide. Currently, the speed brake mode is turned off for our model. 

Three different command signals, which must be converted into excitation signals for the two 
pairs of control surfaces, come from the D1 controller. The following flow chart (Figure 4.3) illustrates 
the methodology used to blend the control surface movements. Basically, this control surface 
management logic assigns priority to pitch axis control over the other two axes. 
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Figure 4.3 Elevon control management logic flow chart 


Since optimization of the control surface allocation is beyond the scope of this study, a very 
simple surface management logic is developed and employed here. However, this allocation is an 
important aspect of the ultimate design of the flight control system. An important question yet to be 
answered is whether the DI design process flows naturally into this optimization scheme, or whether it 
makes optimization more obscure. Figures 4.3 (above) and 4.4 (below) illustrate control surface 
management logic for rudders and elevons. 
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4.2 Sensor Modeling 


As in other aircraft, direct measurement of sideslip angle on the X-38 is highly uncertain. So, 
the sideslip angle is estimated by combining other measurements. The following equation is one way to 
estimate sideslip angle. An integral form of this equation is used in simulations for this study. 


r meas + 


ymeas 


6 meas + 


ymeas 


All other parameters are assumed to be measurable. White noise is added to the outputs to mimic 
sensor noise. 


4.3 Gust Modeling 

Gust effects, modeled as a disturbance, are also considered when evaluating the DI controller. 
Gusts are added to the body component velocity of the vehicle as shown in the following block diagram 
(Figure 4.5). Typical gust inputs are also shown in Figure 4.6. 


Gust Model Input 


Body Component 
Velocities 


Output Velocities 


Figure 4.5 Gust modeling 



Figure 4. 6 Typical gust inputs 
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4.4 Comparison Between MACH Controller and TAMU 
Design 


4.4.1 Control Variable Definition 

Derivatives of the state variables p , q, and r are used as CVs for the X-38 DI example. Instead 
of specifying p, q. or r, however, MACH specifies the dynamics of the derivative of the CVs. The CV 
vector has three components: LCV\ MCV, and NCV. Each of these components controls roll, pitch, 
and yaw moments, respectively. The variables are defined as follows: 


CV = 


LCV 

MCV 

NCV 


(4.7) 


with 


LCV = Ps 


(4.8) 


MCV =q+ 


X 


nz 


nz + - 


Vcois + ^nz) Vcos $ 


-cosycosp.- p s tan$+Kqq- cosy 


■ + 


9 


v v , ; 


(4.9) 


co 


9 


NCV = r s + Kp| 3 - ^cosysinv 


(4.10) 


In these equations, ji represents the bank angle about the velocity vector. 

Since the CVs are no longer states themselves but are combinations of states, a CV definition 
block was added to our Simulink example to “upgrade” to the MACH- generated controller. Although 
differences in roll axis CVs are minor, they can occur in either body-axis roll rates with flight path 
components or in stability axes. And although the definitions in pitch and yaw axes CVs for the MACH 
controller differ from this example, both have a strong angular rate content. Because the first terms on 
the right side of each equation dominate, MCV and NCV essentially become q and r, respectively, in 
our example controller. 

4.4.2 Desired Dynamics Module 

The Desired Dynamics module of the MACH controller is, in form, a PI controller that 
is identical to our Simulink example in Figure 3.10. It has the same magnitude of bandwidth, A#, of 5 
sec ' as well. 


4.4.3 Dynamic Inversion 

The MACH controller starts with the same nonlinear form for the vehicle equations of motion as 
outlined in Section 3.2. 


j&= F(u,x ) 

~f{x) + g{*)u 


(4.11) 
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where x is the vector of state variables and u is the vector of control effectors. Recalling that the 
CVs in MACH are functions of the state variables, 


CV = CV(x), 

where CV is the CV vector defined in Equation (4.7). Then, 


■ dCV „ dCV , dCV . . 

cv= ^r^^r ,{x)+ ^r B{x)u - 


So, the DI control law can be obtained as 

pcv 

dx 


(4.12) 


(4.13) 


g(x) 


- 1 


• <tes i)CV 

CV -^Cf(x) 

dX 


(4.14) 


4.4.4 Control Effector Priority (Surface Management) 

As was previously explained in Section 3.2, the first bracket in the right side of the Dl control 
law must be a square matrix so it can be inverted. This mathematical requirement means that the 
number of control effectors must equal the number of CVs. In the MACH application, the row 

dimension of is 3; therefore, the column dimension of g{x) must be 3 as well to have a square 

dx 

matrix. Thus, three control effectors are required. Since the three CVs in MACH have a strong angular 
rate content and the X-38 lias three moment -producing controls (differential elevons, symmetric 
elevons, and rudders), the vehicle has an appropriate number of control effectors for three of the four 
flight control modes listed in Table 4.3. 


Table 4.3 MACH V201 Flight Control Modes 


Mode 

Flight 

Mode 

Conditions 

#of 

CVs 

CV Feedback 

#of 

Effectors 

Effectors 

1 

ACS Only 

q 'bar < 2 

3 

INS: p. q, r, p 

3 

T , T v , T_ 

2 

Blended 

ACS 

Kj 

A 

A 

3 

INS: p, </, r, p, N : 
NAVDAD: a, V, q h „ r 

5 

T„ T, T, 8,„ 
5, 

3 

Flaps 

Only 

cji,w >20 
and M < 6 

xL 

INS: p, q , r, N- 
NAVDAD: a ,V,q lH „. 

2 

5 n , 6 ( . 

4 

Rudders 
and Flaps 

M >6 

3 

INS: p, q. r, N : , N y 
NAVDAD (or FADS): 

i tfhnr 

3 

8«, 5,, 5,. 
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Note: T x , T y , and T z are* the reaction control thrusters , cm</ <5* 5* and 5, are the resulting 
effective control surface deflections. 


The four flight control modes presented in Table 4.3 fit two different situations (for the MACH 
controller) as follows: 

1 . The number of controlled variables equals the number of available control effectors (modes 
1, 3, and 4). 

2. The number of controlled variables is less than the number of available control effectors 
(mode 2). 

In mode 3, there are only two available control effectors; the rudders are fixed. Consequently, 
the number of controlled variables is greater than the number of available control effectors. To 
overcome this problem (having more unknowns than equations), the number of controlled variables is 
reduced to two by combining roll and yaw signals. Then, the DI control law can be obtained. 

When redundant actuators exist (mode 2), a unique solution for the actuator commands can be 
found by effectively reducing the number of available controls. MACH applies a weighted and biased 
pseudo- inverse solution. Details of this approach are spelled out in the literature 1 '. 

The MACH controller also incorporates an algorithm called sum that deals with all possible 
combinations of the number of unknowns and the number of equations described in this section. The 
sum algorithm is described in detail in Appendix A of Reference 10. 

4.4.5 Least-Squares Aerodynamic Model 

Tire Simulink- based DI control examples presented in the following sections use time- invariant 
coefficients. Unfortunately, aerodynamic coefficients vary over time for all practical applications. All 
coefficients in the DI controller must therefore be updated to account for flight conditions in which the 
vehicle is operating, and how these conditions change the aerodynamic coefficients. However, it is not 
practical to carry all aerodynamic data at every individual point throughout the entire flight envelope. 

For purposes of making appropriate updates, least- squares curve- fitting functions are applied to the 
aerodynamic data a priori. This curve- fitting function is also carried on board the aircraft to generate 
the required coefficient updates. The six aerodynamic coefficients used in VI 32 MACH modeling are 
functions of the following variables: 

Longitudinal: C,- = Kj (a)-^—q + K 2 (cx)8 e + K 3 (cx) , i = D,L,m 

L L 

Lateral: C/ = K 4 (a)p + K 5 (a) — p + K^ (a)— r + K 7 (a)8 a + K 8 (a)8 r , i = Y,l,n 

All Als are a function of cx only, and the Mach number is held constant at 0.6 in the described MACH 
controller. 

When the same set of coefficients is applied to V201 in the TAMU example controller, the 
least-squares curve fitting is modified to be nonlinear both in oc and in the Mach number for longitudinal 
and lateral/directional coefficients. 

Lateral/directional derivatives are fonned using an elevon trim angle setting, 8 C , based on a linear 
relationship of a for longitudinal trim over the Mach range. 
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Longitudinal: 


C k = C k (a,M) + C kq (a,M)^+C k&e ( a,M)b e , k = D,L,m 


Lateral: 


C„ = c tp (a, M,5 e )P + C kp (a, M,5 e )|£ + C„ r (a, M,8 e ) A 


+C/( ga ( a, M, 5 e ) 5 a + C hr (a, M,d e ) 5 r , k = Y, I, n 


with 5 e = m(M)a + b(M) 

4.4.6 Outer Loops 

4. 4. 6. 1 Bank angle outer loop 

Block diagrams of the bank angle outer loop used in the Simulink example and MACH 
controllers are shown in Figure 4.7. 



Figure 4. 7 Comparison of roll angle outer loop structure 


Obviously, both controllers have the same structure. Also, the outer-loop gains, A' 0 and b { j,, are 
both normally set to 1 sec 1 before they are tuned to achieve desired performance. 

The f$ es signal is passed through LCV Command Limit block for the MACH controller. This 

limiter consists of 

1 . an absolute command limit based on JSC specifications, 

2. a command limit that prevents uncontrollable inertial coupling into the pitch axis, 

3. a command limit that prevents uncontrollable inertial coupling into tine yaw axis, and 

4. an absolute command minimum that allows some commands to get through if the inertial 
coupling limits go to zero. 

More details of the Command Limit block in the MACH controller can be found on pages 4 and 5 of 
Reference 10. 


4. 4. 6. 2 Angle-of-attack outer loop 

The alpha outer loops for the Simulink example and the MACH controller also have essentially 
the same structure (Figure 4.8). Both have unity feedback gain (A",, = h„ = 7). Once again, slight 
differences in the two controllers come from definitions in the command variables (MCV ^ q). 
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Figure 4.8 Comparison of angle-of-attack outer loop structure 


Finally, in terms of outer- loop structure, neither controller contains an outer loop that will control 
yaw angle. 

4.4 . 7 Comparison of Aircraft Models 

The following nonlinear model " is used for the MACH controller: 


l xx &~ /« &= la + Ip + lyy r q + l xzP c l ~ l zzP r 
lyy#=rn c a + m c p - / xz p 2 + l^pr- l xx rp+l xz r 2 
-lxz&+ Izz $=n c a + n c p + l xx pq- l xz qr- lyypq 


V&= — l-Dcosfi + Y sin$ + T cos$ cos a ) - gsin y 
m 

#= - [ Dsin p cos p + Yco s p cos p + Ls/'np + T (s//7 p sin a - cos p sin (3 cos a )] 

mVcos y L 

1 V-. geos y 

^b= — — [-Ds/'n p s/np - Ys/np cosp + Lcos p + T (cosp sin a + s/np s/np cosa)J — 


mV 

pcosa + rsina + r Os/np cosp tany + Ytany cos p cosp + L ( tan p + tany sin p ) 

cosp mV L 

. ,-i geos ycos\i tan p 

+7" ( sin a tany sin p + sin a tan p - cosatan ycosp sin P)J — 

&= q - (anP( poos a + r sin*) - (L - + r »«) + Vcosp ' 

o 1 „ _ . „ . gcosys/np 

$= -rcos a + psin<x + (Ds/np + YcosP -Ts//ipcosa) + — 

ml/ Y 

$= Vcosycos x 
Vcosysinx 
&= Vsin y 


with 


37 






'a = ^SbC, 
m a = qSbCm 
n a = qSbC n 


4.4.8 Sensor Processing 

The MACH controller assumes full- state feedback. However, since the air data system for the 
X-38 vehicle does not provide information regarding sideslip angle p, this angle must be estimated. 
Figure 4.9 shows the block diagram of the sideslip estimation block in the MACH controller. 


P. q. r 


(t>, o 


V,ct 


Onboard 

Model 


S at ^ «i 


$model 


1 _ 

S 


Hp 


P 

-► 

Estimation 
Gain 




n v accelerometer 


n v accelerometer model 


Figure 4.9 Sideslip Estimation block diagram (MACH controller) 


4.5 X-38 Mathematical Model 


4.5. 1 Overview and Vehicle Parameters 

X-38 vehicle coefficients were extracted from the output of a modified version of the SES lh . 
The subroutine X35_AERO_DERJVS calculates aerodynamics coefficients periodically along the 
vehicle trajectory at user- specified sampling intervals. Since this subroutine does not provide trimmed 
values of coefficients, for this study a trim routine was not added to the program. The coefficients used 
in the simulation are instead taken from untrimmed flight conditions. 

Table 4.4 shows some of the parameters of the X-38 V201 vehicle used during the SES 
simulation as well as the DI controller simulation that will be described in the next section. 


38 





Table 4.4. Mass Properties and Geometry for the X-38 


Parameter 

Value 

Units 

Weight 

17578 

lbs 

x-axis inertia, l x 

8090 

slugs/ft' 

y-axis inertia, I y 

25900 

slugs/ft' 

z-axis inertia, I z 

29200 

slugs/ft : 

x-z product of inertia, I x/ 

1300 

slugs/ft' 

Reference area, S 

260 

ft 5 

Reference span, b 

27.5 

ft 

Reference cord, c 

10.8 

ft 


4.5.2 X-38 Equations of Motion 

The following nonlinear set of equations is used and integrated in the simulation: 

^=Lp(5 + L p p + L r r +L 6g 5a + L & 8r 

&=N p P + N p p+N r r + N, a 8a + Nt5r 

q N M a a + M q q+M Se de 

(K 

<$= p + (sin <p )(tan G)q + (cos p \tan 6 )r 
#= (cosp)q -(sinp)r 
y&= ( sin p Xsec 6 )q + ( cos p )(sec 9 )r 
= X v V T + X a a - g(cos Y )9 + X^de 


Y Y 
* p +^p + 

y T 


v T 


r + —p + —dr 
Vr V T 


^V T+ ^a 


9(sin y) 


9 + 




-fie 


Vr r V T 


5e 


0 gcosy i , ,x 
)&= ,, L (cosp- 1) 


V: 


&= V T siny 


(4.15) 


4.6 Design Example 1 

This section shows the simulation results obtained using the D1 controller developed in Sections 
3.4, 3.5, 4.1, and 4.2 of this document for the X-38 vehicle. The structure of the simulation, as shown 
in Figure 3.5, is built up in a MATLAB/Simulink environment. A total of ten different simulation cases is 
presented for three different flight conditions. The objective of this simulation exercise is to illustrate the 
characteristics of the DI controller as applied to the X-38 lifting-body reentry vehicle with its rather 
laige flight envelope. 
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4. 6. 1 Flight Conditions 

Three flight conditions along the nominal trajectory of the X- 3 8 are selected and used as initial 
conditions of the simulations to evaluate the performance of the designed controller. Those flight 
conditions are selected to represent three different regimes of the flight envelope - supersonic, transonic, 
and subsonic. Table 4.5 summarizes the fight conditions evaluated. Since the SES does not provide 
trimmed values of the vehicle aerodynamic coefficients, pitch, roll, and yaw coefficients are nonzero 
values. 


Table 4.5 Summary of Evaluated Flight Conditions 



Ho(kft) 

M„ 0 


0C„(°) 

M°) 



HE51 

r () (°/sec) 

m 

76.2 

2.38 

276.9 

26.9 

-5.6xl0" 5 



0.13 

0.34 

B 

46.4 

1.05 

221.9 

16.3 

-1.3x1 O' 2 

-1.68 

-5.2x1 0 3 

-5.8xl0' 2 

-0.12 

H 


ESI 

267.9 

11.9 

1.6x1 0‘ 2 

0.0 

1.0x1 O’ 2 

-0.43 

5.9x1 0‘ 2 


4.6.2 Simulation Run Matrix 

Ten runs were made to demonstrate the characteristics of the DI controller. This set of 
simulations is summarized in Table 4.6. The test matrix consists of three subsets of test objectives. The 
first set. Runs 1 through 4, shows the vehicle responses during typical maneuvers performed in each 
portion of the flight envelope using the nominal set of aerodynamics coefficients. The second set, Runs 
5 through 9, shows how uncertainties in aerodynamic coefficients affect the performance of the DI 
controller. Finally, Run 10 shows the effectiveness of the controller in coping with a side gust (i.e., 
acting essentially in the xy-plane of the vehicle), an outside disturbance for this mathematical model. 


Table 4.6 Simulation Run Matrix 


Run No. 

Flight Cond. 

input 

Remarks 

I 

A 

-5° a- Step, 25° <f>- Singlet 

Nominal Case 

2 

B 

-0.05°/sec a-Ramp,7.5°0-Singlet 

Nominal Case 

3 

C 

7.5° (J)-Singlet 

Nominal Case 

4 

A 

7.5° 0- Singlet 

Nominal Case with FQ, = 0.4 

5 

A 

10° a-Singlet, 10° 0-Singlet 

Nominal Case 

6 

A 

10° a -Singlet, 10° 0- Singlet 

30% Aero Uncertainty* 

7 

A 


50% Aero Uncertainty* 

8 

A 

ft B!S BM 

60% Aero Uncertainty* 

9 

A 

10° a-Singlet, 10° 0-Singlet 

60% Aero Uncertainty*, K <i) =0.4 

10 

C 

Side Force Gust (50 ft/sec max) 

Side Force Gust 


*See Table 4.7 for complete definition. 
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Table 4 . 7 Aerodynamic Uncertainty Matrix 


























































Run 1 (Figure 4. 10) is the nominal case for supersonic flight conditions. A rather aggressive 
step input of -5° in a-command and a 25° singlet input in ^-command were used as command 
excitations. Both inputs were engaged at 1 second, and the duration of the singlet was set for 9 
seconds. Despite the aggressive commands, the vehicle behaved very well both longitudinally and 
laterally/directionally, with little overshoot about any axis. Also, no actuator saturation occurred with 
either rudders or elevons throughout the simulation. 



0 10 20 30 0 10 20 30 



Left Elevon (deg)0 20 30 Right Sevon (de1$ 2 0 30 

Time (sec) Time (sec) 


Figure 4.11 Simulation Run 2 , transonic flight 


Run 2 (Figure 4.1 1 ) demonstrates nominal transonic performance. This time a ramp of only - 
0.05 deg/sec for a and a 7.5° singlet for (j) were used as command inputs. Despite the fact that the 
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controller was unchanged from the previous case, the vehicle again responded very well and no 
actuators saturated. 


A' 




Figure 4. 12 Simulation Run 3, subsonic flight, original unity outer loop gain 


Run 3 (Figure 4. 12) illustrates performance of the DI controller in subsonic flight, a 
was commanded to maintain its initial value while <j) was commanded to follow a 7.5° singlet. The 
controller functioned well until about 10 seconds into the simulation, when elevon saturation occurred 
and the performance degraded rapidly especially in the roll axis. The roll axis degraded more than the 
pitch axis because, as explained previously, the control surface management logic puts priority on the 
pitch axis over the roll axis. There is little augmented damping in roll during the elevon saturation period. 
To overcome this deficiency, the gain of the bank angle outer loop was reduced from 1 to 0.4, and the 
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simulation was repeated using the same inputs. Output from this modified controller is shown in Figure 
4. 13 (Run 4). Although the elevons still saturate, the overall performance is improved over the unity 
feedback gain case. 
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Figure 4.13 Simulation Run 4, subsonic flight, outer loop gain = 0.4 


One point of emphasis: Although D1 advocates usually assert no gain scheduling is needed 
because of DFs “global linearization” structure, in reality gain adjustments are still required whenever 
commands are large enough to saturate actuators. DI performs poorly during actuator saturation 
because only the aircraft dynamics are inverted and actuator dynamics are simply left out of the 
controller formulation (at least in this implementation). In other words, some optimization scheme that 
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recognizes actuator constraints is needed to overcome saturation problems. This drawback has more to 
do with the control surface allocation scheme than it does with the D1 design itself. However, this result 
clearly suggests that the control allocation issue needs to be attacked by treating actuators as nonlinear 
elements and including them in the DI process. This effort will be handled in a future study. 

4.6.4 Uncertainties in Aerodynamic Coefficients 

The next issue to be discussed is uncertainty in aerodynamic coefficients. Out of the many 
parameters used in the simulation, 1 3 major aerodynamic coefficients were selected; they are presented 
in Table 4.5. The uncertainties are represented as algebraic additions to each parameter, and the 
magnitude of the uncertainties is set to three different percentages: 30, 50, and 60. The signs of these 
uncertainties, either an addition or a subtraction to/ffom the nominal values, are randomly chosen and 
are also shown in Table 4.7. 

This set of simulations consists of five runs, one run for each nominal case as well as 30% and 
50% uncertainty cases and two runs for the 60% case. The same set of flight conditions (A in Table 
4.6) is used for all five runs. The same set of command inputs, 10° singlet commands for both a. and <>, 
is also applied for each run. 

The controller used for Runs 5 thorough 8 is the same controller that is used for Runs 1 
thorough 3. The controller used for Run 4 has the same structure, but the gain of its outer 4)- loop is 
reduced from 1 to 0.4. This modified controller also is used for Run 9. 

Run 5 (Figure 4. 14) is the nominal case. The vehicle responded well to the commanded inputs 
and no actuator saturation took place. 

Run 6 (Figure 4.15) is the 30% uncertainty case. Despite the right elevon saturation at 
2 seconds into the simulation, the vehicle still behaved well. 

Run 7 (Figure 4. 1 6) is the 50% uncertainty case. The divergence from the nominal case is 
obvious. Overshoots are large, especially during the period when the elevon saturation occurs. The 
bank angle excursion is also large because of pitch axis priority in the control surface management 
scheme. The lack of roll control authority under these conditions is underscored in this simulation. 

Next, uncertainty is increased to 60% in Run 8 (Figure 4.17). The vehicle obviously diverges 
from the nominal trajectory, probably catastrophically. 

From the 50% uncertainty case results, as well as by looking at the right elevon time history in 
this case, it is clear that elevon saturation, coupled with the priority for pitch axis control, causes loss of 
control. A quick (though not necessarily optimal) solution to this problem is, as in the previous case, 
reduction in outer (j)-loop gain. This approach is examined in the next simulation run (Figure 4. 1 8). 
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Figure 4.15 Simulation Run 6 , supersonic flight (M„ = 2.38), 30% uncertainties 


7 












40 


2 






Left Elivon (degt)0 20 30 

Time (sec) 




Ri ght f§levon (de^ 2 

Time (sec) 


Figure 4.16 Simulation Run 7, supersonic flight (M„ = 2.38), 50% uncertainties 


By using the controller from Run 4 with the outer 0-loop gain reduced from 1 .0 to 0.4, 
simulation Run 8 was repeated as Run 9. With this change in gain, the vehicle readily recovers from the 
upset, which suggests that gain modifications are effective in alleviating actuator saturation. Certainly the 
complexity of the controller increases as the number of gain alterations increases. There is, therefore, a 
tradeoff that designers have to consider between controller complexity and performance of the 
controller for the best controller design. 


48 













Left EI<fVon (deg)IO 20 30 

Time (sec) 


5 
0 
-5 

Beta^ 

0 10 20 30 




iffi 


m 



IB 


Ri ght E&von (deg^ 20 

Time (sec) 


30 


Figure 4.17 Simulation Run 8, supersonic flight (M„ = 2.38), 60% uncertainties 
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Figure 4.18 Simulation Run 9, supersonic flight (M„ = 2.38), 60% uncertainties, 

outer O-loop gain = 0.4 


4.6.5 External Disturbances Effect: Side Gust 

The final simulation example is an external disturbance case. We chose a side gust as a 
representative case of such disturbances and their effect on D1 controllers. The gust input shown in the 
upper left comer in Figure 4. 1 9, Run 1 0, was added to the body y-axis velocity of the vehicle. The gust 
reaches up to 50 ft/sec, and the overall duration of the gust is 4 seconds. The flight condition used is the 
subsonic case (M„ = 0.63), Flight Condition C, shown in the Table 4.6. 
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The side gust disturbance input produces a maximum sideslip angle of about 1°. Also, bank 
angle is excited up to about 1 2° because of yaw- roll coupling. Bank angle did damp out rapidly (after 
approximately 1 cycle), and the sideslip angle became essentially negligible roughly 10 seconds after the 
gust was introduced. Overall, this DI controller worked well for rejecting external disturbances. 









Figure 4.19 Simulation Run 10, subsonic flight (M„ — 0.63), 
external disturbance: side gust 
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4.7 Design Example 2 


4. 7 . 1 Introduction 

Even in instances where the control input matrix g(x) shown in Equation (3.4) is invertible, 
problems will arise if the matrix is very small in magnitude, indicating reduced control effectiveness on 
the state dynamics. Should this occur, the control becomes unbounded, causing actuator saturation. 

For example, a control ineffectiveness problem is typically encountered in the inversion of angle-of- 
attack dynamics because the vehicle’s longitudinal control surface has little effect on the angle-of-attack 
rate. A 2- time scale method has been developed and applied in previous research to bypass this 
problem. 17 18 In our method, the control surface is used to generate the pitch rate dynamics, <$, 
directly. The resulting pitch rate is then used to control <&. Figure 4.20 illustrates this approach in which 
two DI loops are present: a fast inner- loop inversion for rotational rate variables, such as q, and a slow 
outer- loop inversion for rotational variables, such as a. 



Equation (4.16) shows the pitch rate command generated by the slow outer- loop inversion, and 
Equation (4.17) shows the elevon command generated by the fast inner- loop inversion, required to 
produce the desired angle-of-attack response. It should be noted that both control laws require angle- 
of-attack and pitch rate feedback, which is full -state feedback for the short period approximation. 


Qcind — 


' Z - Y VK 


i+-s 


v 


8 c , m> .=M- | (C.-M„a-M q q) 


(4.16) 

(4.17) 


The control laws shown above are used to invert the inherent pitch rate dynamics and the angle-of- 
attack dynamics, respectively. The proposed forms of desired dynamics developed in Section 3.6 are 
applied separately to each of the short-period state dynamics before they are evaluated in terms of 

• Time domain performance, 

• Stability robustness and performance robustness, 

• Effect on motion sickness, 

• Quadratic cost function, and 

• Passenger ride comfort index. 
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4.7.2 Time Domain Design Requirements 

The time domain performance specifications used to evaluate vehicle response are short-period 
damping ratio and natural frequency. The requirements for the X-38 vehicle are selected from Mil- 
STD- 1 797 A 14 to satisfy Level 1 flying qualities for a Class II vehicle during a Category B flight phase. 
These Level 1 standards bound the short-period damping ratio and natural frequency as follows: 

0.30 <£< 2.0 (4.18) 

0.7s" 1 <(0 n < 1.3s" 1 (4.19) 

The time domain response requirements for a step input are shown pictorially in Figure 4.21. These 
requirements satisfy Level 1 flying qualities for the vehicle. 



0 T i T 2 T 3 T 4 T 5 T 6 



Figure 4.21 Time domain performance specifications 
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Two numeric indices are also used to evaluate the responses of each controller. A quadratic 
cost function, J, (shown in Equation (4.20)) is used to evaluate the magnitude of each controller 
response in terms of the four longitudinal states and the elevon input. 

J=x T Qx+u'Ru (4.20) 

It was desired to weight angle -of- attack more heavily than pitch rate because angle -of- attack is the 
variable being controlled. Since the behavior of the phugoid mode is less important, the velocity and 
pitch attitude states are weighted less than the short period states. Therefore, the weighting matrices 
used are: Q = diag(.l 10 1 ,l)andR=l. A passenger ride comfort index 19 , shown in Equation 
(4.2 1 ), was also used to evaluate the responses. 

C = 2. 1 + 1 7.2a, (4.21) 

The RMS vertical acceleration, a z ,is a, = U 0 fife-U 0 cos(a o )q+gsin(6 o )0 in which vertical 

acceleration subscripts denote trim values. Lower comfort indices indicate increased passenger ride 
comfort. 


4.7.3 Controller Design 

The proportional, PI, and flying qualities dynamics desired are selected to satisfy Equations 
(4.18) and (4.19). The ride qualities dynamics were selected to satisfy the desired CAP and damping 
ratio values shown in Figure 3.13. Table 4.8 summarizes the desired dynamics selected for these 
controllers. Each set of dynamics acts on the error between CV command and its feedback term. The 
robust outer loop was designed via pole placement through the use of MATLAB’s place command. 
This outer loop consists of full- state feedback that is operated on by a matrix of gain values. The poles 
were placed according to the desired closed- loop pole locations. 


Table 4.8 Desired Dynamics Selection 


Desired Dynamics 

Angle-of-Attack Case 

Slow inversion 

Fast Inversion 

Proportional 

0.8 

1.3 

Proportional 

Integral 

( 1 ^ 

6 — CV -CV 

) 

+— (CV_-CV) 
s 

(\ \ 

— C\! -CV 

\ z y 

+-(CV_-CV) 

s 

Flying Qualities 

1.4(s + 0.8) 

1.5(s + 0.8) 

s' + 2.24s + 1.96 

s' + 2.24s + 1 .96 

Ride Qualities 

1.96 

S + 3 

7 

s + 4 
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4. 7.4 Time Domain Analysis 

Time histories for the inverted angle-of-attack d>namics are shown in Figure 4.22. The PI 
dynamics respond very fast, causing a large actuator rate initially. Although the flying qualities dynamics 
slightly violate the angle-of-attack time domain constraints, the other forms of dynamics stay within the 
time-domain specifications. It should be noted that all responses satisfy the actuator position and rate 
constraints. Of the various desired dynamics, the proportional form results in the lowest cost and lowest 
comfort index for an angle-of-attack DI controller. The resulting costs and comfort indices of these 
responses are shown in Table 4.9. 
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Figure 4.22 Time histories for the inverted a dynamics 


Table 4. 9 Cost and Passenger Comfort Index 


Desired Dynamics 

Cost 

Comfort Index 

Proportional 

130.4 

2.56 

Proportional Integral 

1 72.0 

2.60 

Flying Qualities 

142.4 

2.58 

Ride Qualities 

139.7 

2.61 
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4.7.5 Frequency Domain Analysis 

Robustness is a key element in the reentry vehicle flight control system design because of the 
broad flight envelope in which these vehicles must operate. The robustness technique that is used to 
analyze controllers developed in this work is adapted from current research in robustness methods for 
the X-38 vehicle 20 . A sigma- Bode plot of the loop gain singular values is used to evaluate robustness 
over a range of input frequencies. Further explanation of the sigma-Bode plot can be found in Section 
4.8.2. The performance and robustness criteria used to evaluate the controllers in this section are 

• Zero steady- state error. 

• Attenuation of low-frequency disturbances by a factor of 0. 1 . 

• Linear model accuracy to within 10% of actual plant for frequencies up to 2 rad/sec and growth 
without bound at 20 dB/decade thereafter. 

To satisfy these requirements, illustrated in Figure 4.23, the singular values of the loop gain must lie 
outside of the performance requirement and stability requirement areas for all frequencies. 
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Figure 4.23 Robustness constraints 
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An additional requirement on signal attenuation can be added to the singular value plot of Figure 
4.23. For this example, attenuate the signals between 0.6 rad/sec to 1 .6 rad/sec in order to alleviate the 
passenger motion sickness that can occur at this frequency range 10 Since the singular values of a matrix 
represent the relative size of a mode, attenuate the singular value response within this frequency range. 
The singular value response within this range should be continuously decreasing. Any amplification 
within it would increase the passengers’ motion sickness. 

The sigma-Bode plots for these cases are shown in Figure 4.24. Both the proportional and the 
PI desired dynamics violate the stability robustness requirement above 20 rad/sec. These dynamics also 
violate the attenuation requirement toward the higher end of the frequency range. Figure 4.24 clearly 
indicates that flying qualities compensation and ride quality compensation are sufficient for a robust 
angle-of-attack DI controller because they satisfy stability robustness, performance robustness, and the 
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motion sickness constraint. Table 4. 1 0 summarizes the results of each set of desired dynamics in terms 
of satisfying time domain and frequency domain requirements. 



Figure 4.24 Sigma-Bode of closed-loop system 


Table 4.10 Summary of Compliance with Design Specifications 


Desired Dynamics 

Step 

Response 

Control 

Responses 

Robustness 

Constraints 

Motion 

Sickness 

Proportional 

3 

3 

5 

5 

PI 

3 

3 

5 

5 

Flying Qualities 

5 

3 

3 

3 

Ride Qualities 

3 

3 

3 

3 


4.8 Design Example 3 

4.8.1 Introduction 

A systematic way to synthesize and analyze the robustness of a Dl-based controller 
in a straightforward manner must be found so that engineering intuition can be easily applied throughout 


57 







the design process. The approach presented in this study uses LQG techniques 21 to synthesize the 
outer- loop controller. This approach is a relatively easy-to-use design method for multivariable control 
design. It is a loop- shaping tool in the frequency domain that gives robustness to the system 22 ' 25 . 
Desired dynamics are given by a dynamic compensator that shapes the loop. Selected dynamics are 
based on performance and stability robustness requirements. These requirements are straightforwardly 
formulated during synthesis of the controller as frequency-dependent singular value bounds (Figure 
4.25). 



Figure 4.25 X-38 lateral-directional control system 


4.8.2 Design Requirements 

The controller design procedure is iterative and centers on designing a controller that will satisfy 
a set of specifications. Specifications are usually stated in both time and frequency domains. A time 
domain requirement for the X-38 lateral-directional system is defined by the time histories of a 10° bank 
angle response after a step input 26 . The boundaries of this requirement are depicted in Figure 4.21 . 

Two different Mach- number- dependent requirements are specified in the time domain. Clearly, 
requirements for the subsonic regime are much tighter than for supersonic flight conditions. According 
to the requirements 26 , no sideslip angle constraint is spelled out for lateral-directional vehicle control. 
Nevertheless, there should be as small a sideslip perturbation as possible to maintain a coordinated turn 
at all times. The D1 inner- loop controller cancels the existing system dynamics and replaces them with 
designer- specified responses. Since desired outputs are often decoupled about each axis, off-diagonal 
coupling effects are typically minimal for a good DI controller. 

Next, frequency domain specifications are meant to ensure performance and stability 
robustness. Both of these requirements are expressed using singular values. Singular value is a suitable 
choice to express the magnitude of matrix functions because it generalizes known SISO statements and 
constraints of the design problem to MIMO cases. In general, a singular value is thought of as the Bode 
magnitude plot for an SISO case, but singular values extend the concept to a MIMO system as well. 
Therefore, MIMO design is carried out using classical control design concepts that align with 
engineering intuition. 

To achieve an adequate response, the gain in the low-frequency region must first be high enough 
to give a quick response to the input while the slope of the singular values must be steeper than -20 
dB/decade to reduce the steady-state error. Further, we assume it is desirable to have at least 0. 1 
rad/sec of crossover frequency to obtain a good closed- loop transient response. These requirements 
are integrated to form the singular value bounds in the low-frequency region. This low-frequency 
“trapezoid” is sketched in Figure 4.26. 
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Frequency (rad/sec) 


Figure 4.26 Frequency domain requirements 


At high frequencies, singular values are bounded by the unmodeled dynamics associated with 
high-frequency models such as flexible and vibrational models. These high-frequency models are often 
neglected when the plant model is being developed; the vehicle is often treated as a rigid body. As a 
result, there is a difference between the assumed mathematical plant, gDC, and the actual plant, (7 'DC. 

This difference is defined and described as model uncertainty 27 . Here, an unmodeled dynamics model, 
suggested by Stevens and Lewis 27 , is used. This uncertainty model assumes the rigid body model is 
accurate to within 10% and up to a frequency of 2 rad/sec, after which the uncertainty grows at a rate 
of 20 dB/decade. This uncertainty model is expressed in the following transfer function and is used to 
model the uncertainties of the X-38 vehicle at high frequencies: 

,„bg^ , 4 . 22 , 

We assume ///DC to be bounded with uncertainty in the X-38 vehicle transfer function. Multiplicative 
uncertainty is expressed in terms of an assumed plant model and the actual plant by 

G'B»Q[/ + (4.23) 

where the unknown discrepancy satisfies a known bound. 
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crO/^o^lmD( 


(4.24) 


Under this assumption, for stability robustness with modeling errors, the loop gain referred to the output 
should satisfy: 


g(&K 



(4.25) 


when l/«70Q:<l. Here, K represents the compensator. This uncertainty bound, which is obtained 

from the unmodeled dynamics, is applied in the high-frequency region. Finally, the complete frequency 
domain bounds are shown in Figure 4.26 (previous page). 


4.8.3 Lateral-Directional Dynamic Inversion Controller 

Lateral- directional D1 control equations were developed previously in Section 3.2.2. The 
following form, which is provided in Equation (3.12), is used for this example: 


da 

and 

A* 

-1 

<{ 

"A 

dcs 

'K A 

dr 


N Sr _ 


_& 




N a 


P 

r 

P 


(4.26) 


4.8.4 Dynamic Inversion Inner-Loop Controller 

The lateral-directional D1 controller, together with the X-38 model and control surface blocks, 
forms the Df inner- loop augmented system. This augmented inner loop is 

ifyj — \ n X U + ByyU^ , (4.27) 

where the inner- loop state vector, x n , is 

x lL = [(3 p r (}) do8c da && <5r] 7 (4.28) 

and the inner- loop control vector, u„ , is 


u 


IL ~ 



(4.29) 


The state matrix, A ;/ , and the control distribution matrix, for transonic flight (Table 4.5, Case B) 
are expressed numerically as follows: 
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We assume measurements are limited to bank angle, 6, and to sideslip angle, (3. Therefore, the 
following output equation results are: 


A 


'0 

0 

0 1 

0 

0 

0 

0“ 

A 


i 

0 

0 0 

0 

0 

0 

0 


(4.30) 


The block diagram representation of this DI inner-loop is shown in Figure 4.27. In addition, the singular 
values of the inner- loop versus frequency are plotted in Figure 4.28. 


pt* 



Figure 4.27 Dynamic Inversion Control Inner-Loop block diagram 
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Frequency (rad/sec) 

Figure 4.28 Singular values of the dynamic inversion inner-loop system 


4.8.5 Augmented System 

Since DI alone does not achieve the desired specifications, stability robustness requirements at 
high frequencies are not met. Integrators are added to each control channel to correct this deficiency. 
The X-38 plant, actuators, and DI controller are augmented to form the following system: 


= A + B 




(4.31) 


where 


X ( ,„„ =[)3 P r </> 8& 5a 8& Sr e 0 


(4.32) 


and 


u 


au- i! 



(4.33) 
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In Equation (4.32), e 0 and e„ are augmented states because of the addition of the integrators in the 

bank angle and the sideslip channels, respectively. The augmented matrices, A, mi , and are 
obtained by perfonning the following manipulations: 


A = 


X 


and B = 

o' 

0 

0 

any 

I : 


(4.34) 


The corresponding singular value plot for this augmented system is shown in Figure 4.29. It suggests 
that the nominal design has now been reshaped to meet the specified frequency domain requirements. 
But we still have not tailored the outer loop for robustness, and we have not dealt with the LQR-LQG 
observer issue. 



Figure 4.29 Augmented system singular values 


4.8.6 Observer Design 

The objective of this step is to create a fast dynamics observer that can be used with a regulator 
to fonn an LQG controller that will satisfy both performance and stability robustness requirements. The 
form of the observer is a Kalman filter designed for the augmented system of the previous section. The 
following weighting matrices were selected, by trial and error, so that the singular values of the resulting 
Kalman filter open- loop gain, CF (s)L , satisfy the singular value frequency domain requirements: 
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Q = diag( 1, 1, 1, 1, 1, 1, 1, 1, lxlO 8 , lxlO 8 ) and R = 0.0051, 


where F (s) = (.si - A) 1 . 

The resulting Kalman filter gain is given by 

-5. If- 1 2.6f+l 8.7f-l 1.6f + 1 1.4f + 1 i.9 2.6 -8.3f-l 1.4f + l 5Ae-\ 

L = 

1.6f + 1 - 1 ,6f + 1 -2.5f + 1 - 5. If- 1 1 ,4f + 3 6.9f + l -3.2e+2 1.8 5.1f-l -1.4e + l 

and the corresponding Kalman filter poles are 


.y = - 1 8.4 ± 1 8.4, - 1 8.4 ± 1 8.5, - 14. 1 ± 0.72,-0.86 ± 0.50, - 0.87 ± 0.49 . 


4.8 . 7 Regulator Design 

A regulator is designed next, assuming full- state feedback. The resulting regulator is combined 
with the Kalman filter from the previous step to form an LQG controller for the system. The regulator is 
based on LQR methodology, and the following weighting matrices are selected for the LQR gain 
calculations so that the corresponding LQG singular value plot (Figure 4.30) satisfies the singular- value 
frequency domain requirements: 

Q =diag(l 1,1,U,LU, lx 10 s , 1x10 s ) and R = 1, 



Frequency (rad/sec) 

Figure 4.30 Singular values of the LQG regulator 
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The resulting regulator gain is as follows: 


-l.lt?— 1 8.2 1.2 3.0i— 3 3.8f-2 1.4 

-4. le- 2 1.2 3.5t?+ 1 -l.k-1 -4.3t>-3 -2.lt— 1 


3.9c -2 1.4 1.0c +4 

-4.2c -2 -1.5 1.7c- 4 


1.7c -4 
1.0c + 4 


(4.35) 


For the LQG singular value plot, the complete dynamics are given by the following augmented system: 



A 

-BK Tx" 

+ 

0 

0 

A-LC-BK_[x_ 


-L 


y = [c 



X 

X 


(4.36) 


with transfer function CF (s)BKF , (s)L . The singular values are plotted in Figure 4.30. The singular 
values of the resulting system clearly meet frequency domain specifications at all frequencies of interest. 

4.8.8 Time Domain Analysis 

Once the controller is designed to meet the frequency domain specification, its perfonnance 
must be tested in the time domain against the time-domain specification (Figure 4.2 1 ). Figure 4.3 1 
shows the response to a 10° bank angle step input to the system as defined in Equation (4.36). The 
bank angle response is within the design envelope while the sideslip angle is negligible. Also, both 
control surfaces are relatively inactive, and both displacement and rate are within the limits for each 
surface. The controller design, therefore, satisfies time-domain specifications at this stage. 

4.8.9 Gain Scheduling Issues 

As previously discussed, the main advantage over classical methods for the Dl design 
methodology is the little need associated with the Dl design methodology to schedule gains. In order to 
verify this claim, the designed controller, which is tuned at the transonic flight condition (Table 4.5, Case 
B), is applied to other flight conditions (Table 4.5, Cases A and C) without modification. The other two 
flight conditions, subsonic and supersonic, illustrate this analysis. The same 1 0° bank angle step inputs 
are applied to all three cases, and the resulting responses are presented in Figure 4.32 for the subsonic 
and transonic flight condition cases. 
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time (sec) time (sec) 

Figure 4.32 10° bank angle step response for different flight conditions 


In the transonic and subsonic cases, the designed controller is able to stabilize the system. However, it 
fails to stabilize the system for the supersonic flight condition. In transonic and subsonic flight, control 
surface activities are well below the limits in rates and displacement. Though the objective of the DI 
controller is to produce a desired response at all flight conditions, no constraints are imposed on control 
surface activities in the DI control equation to achieve this. Control surface deflections and rates are 
based solely on the control distribution matrix. In order to avoid actuator saturation, the DI controller 
must command no more deflection or rate than the system hardware can provide or the system 
becomes nonlinear and the linear analysis may break down. But, these constraints overly restrict the 
available control power in some cases. Since a controller selected with this process may not produce 
the “best” performance at all flight conditions, the “nominal” (design) condition must be selected 
carefully to achieve good performance in a wide range of flight conditions. Moreover, the specification 
will change for different flight conditions. For example, the time-domain specification changes drastically 
near Mach = 1 . Changing requirements for the entire flight envelope must be considered during the 
design process. The proposed approach allows an engineer to address this issue rather intuitively. 
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5 Robustness Analysis 


5.1 p -Analysis Applied to the X-38 

5.1.1 Introduction 

In the introductory sections of this document, singular- value decomposition is mentioned as the 
most common approach to adding robustness to a DI controller design. But it should be again 
emphasized that p- analysis is by no means the only way to tackle robustness issues with this flight 
control design. Dang Vu 28 suggests that combined techniques using linear quadratic design, quantitative 
feedback theory, Lyapunov synthesis, adaptive control, and differential games have all been initiated 
together with p- analysis. 

5.1.2 Robustness Example: Application to the X-38 Lateral- 
Directional Aircraft Equations of Motion 

The linear fractional transformation (LFT) structure - spelled out in detail in Section 6.2 of this 
document - is now applied to the lateral directional aircraft dynamics for the X-38. First, the lateral- 
directional aircraft equations of motion are represented as follows: 
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(5.1) 


Next, we assume that - except for L p , L&, N r , and A'V - all parameters are certain values. We also 
know the variations (boundaries of uncertainty) for these parameters, too; and we will write them in the 
following form: 


L-<Ln<Ll 


L h * L 6 a * 4 


N~ <N r <N P 
Ns r <N 5r <Nl 


(5.2) 


These bounded uncertainties are now integrated into the aircraft equations of motion using the LFT 
form. Compare the decoupled roll axis equation 

&= L pP + L d a 5 a 

y = p 
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with the general expression in Equation (6.8). Then, the following correspondences are clearly found: 


x = p 
u = b a 

d = Lp 


(5.4) 


Now, ensure the state space representation of the roll axis eqation with bounded uncertainties in L p and 
L&, is written by substituting the preceding relationship into Equation (6.13). 
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(5.5) 



Similarly, yaw axis aircraft equation 


/&= N r r + A/§^5 f 
y = r 


is put into state space fonn with uncertainties in N, and /V&- 
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(5.6) 


(5.7) 
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Integrating both roll and yaw axis equations with Equation (5.1 ) yields the following LFT form for the 
lateral-directional equations: 




0 O' 

1 °\vi 
0 0 v 2 

0 1 V * 


0 0 
0 0 
0 0 
0 0 


(5.8) 


The input/output description of the plant and the uncertainty block are shown in Figures 5.1 and 
5.2. 



Figure 5.1 Plant input/output 
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Figure 5.2 Uncertainty block 


These two blocks are interconnected to form a Parametric Uncertainty block, as shown in Figure 5.3. 
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Figure 5.3 Aircraft plant with parametric uncertainty 


Parametric uncertainties are not the only type of uncertainties to be considered in this example. 
As discussed in Section 6.2.2. 1 , unmodeled dynamics or uncertainty at the input is another important 
type of uncertainty to be examined. When this is applied to our case study, the following weighting 
function, W„„ is used. Its place in the lateral- directional block diagram is shown in Ftture 5.4. 
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(5.9) 
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Figure 5.4 Unmodeled lateral-directional aircraft dynamics 



Frequency (rad/sec) 

Figure 5.5 Uncertainty weighting function 


This particular uncertainty weight, W„„ indicates that there is potentially 1% modeling error at 
low frequency and that the uncertainty in the model grows up to 100% at high frequency. The 
uncertainty weight is diagonal in form with equal diagonal elements. The perturbation mode is a circle, 
or a sphere, around the nominal plant. 

A third type of uncertainty, discussed in Section 6.22.2, is uncertainty at the output (uncertainty 
in the measurements used) in the feedback loop. For our lateral- directional example, Figure 5.6 shows 
a block diagram representing the treatment of this type of uncertainty. Since it is unstructured, it is a 
function of input frequency - as was the uncertainty at the input. 


fimeas 4 " 
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Figure 5.6 Unstructured uncertainty at the plant input due to output uncertainty 
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Frequency (rad/sec) 

Figure 5 . 7 Unstructured output uncertainty weight 


This output uncertainty weighing function implies in p and r a low- frequency measurement error of 
0.003 rad/sec and a high-frequency measurement error of 0.015 rad/sec. The model of measured value 
of /;, denoted p mim , is given by 


Pmeas ~ P + ^n p Hp > 


(5.10) 


where T| /; is an arbitrary signal with 


< 5 "> 

Any type of controller could be used in conjunction with our DI controller for ji-analysis. Much 
of the literature favors an //.. controller that can provide the robustness sought through p- analysis - 
largely because p- synthesis. which is an extended and more complex form of p-analysis, requires an IL 
controller combined with p-analysis to synthesize an optimized controller to achieve stability 
performance. In our example, an //. controller is designed using the MATLAB p- Analysis and 
Synthesis toolbox' 1 . 

The objective of the !L, controller is to make the pitch and yaw rate of a vehicle closely follow 
the commanded values of these two parameters (i.e., the goal is to minimize the errors e,, and e,,). 
Frequency-dependent weights are connected into the structure as shown in Figure 5.8. 



► 6p 


► e r 


Figure 5.8 Performance Weighting block diagram 
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The overall performance error vector is given by 


e perf 


e p 


S + 10 
S+ 1 

0 

Pcmd Pactuaf 

_ e r _ 


0 

s + 10 

_ r cmd ~ r actual 



s + 1 _ 



(5.12) 


The shape of the weighting function (Figure 5.9) is chosen so that the controller provides 
performance in the low- to mid- frequency range. Error weights on the roll and yaw rates indicate a 
tolerance of 0. 1 rad/sec at low frequency and 1 rad/sec at high frequency. 
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Figure 5.9 Performance weighting as a function of frequency 



Limits on the actuator deflection magnitude and rates are also included in this example through 
W ac , (actuator weight) shown in Figure 5. 1 0. This mathematical constraint is not a physical “limit” but is 
treated as a constant weight matrix to produce the “error at the actuator” e,„ t . 



f e acf 


Figure 5.10 Control Surface Actuator Weights block diagram 


This error vector is defined as 
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Next, an IL controller is designed to minimize the H„ norm of the nominal closed- loop transfer 
function from the disturbances to the errors. This design is actually a sub-optimal / /. controller because 
the controller is formulated after assuming that there is no model uncertainty. The D-K iteration process 
(part of the singular value decomposition process) must be done to approach a true optimal / / - 
controller. The input/output relation of the postulated sub-optimal H,. controller is shown in Figure 
5.12. 


Pcmd 
? cmd 

Pmea 
Pmeas 
r mea 


Figure 5.12 controller input/output 
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Although the designed system is a simple system, the order of the controller is 14. High order is 
a disadvantage of M. controllers because system complexity increases as the order of the controller 
goes up. 


“True” Airplane 



Subsystems developed in the previous sections are now integrated to form an interconnection 
structure for p- analysis. A block diagram of the interconnection structure used for the analysis is given 
in Figure 5.13. The main objective in creating an interconnection structure is to transform all of the 
subsystem LFTs to a single (large) LFT that will separate the unknown parts from the known parts of 
the system. The interconnection structure can be formed relatively straightforwardly since standard linear 
operations - such as cascade connections, parallel connections, feedback connections, inversion, and 
frequency response - retain the LFT form. In other words, interconnections of LFTs are still LFTs. 
Since interconnection of the subsystem LFTs is a straightforward, albeit rather tedious, task, using 
commercially available software (e.g., the MATLAB p- Synthesis Toolbox) is recommended. 

Now, we are ready to apply the seven steps outlined in Section 6.2.3 to complete our 
p-analysis and explore the stability robustness for our lateral-directional example. Applying the software 
tools mentioned above to the system shown in Figure 5.13, we complete the p-analysis process. The 
flight condition simulated is the transonic portion of the flight envelope at M, M) = 1 .05, q = 222 lb/ft 2 , 
and a = 1 6.3°. Uncertainty levels vaiy equally among four parameters, L p , Lg , N r and A/g , from 

10% to 1 00%. Figure 5.14 summarizes the results obtained from this p-analysis. 
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Figure 5.14 Parametric uncertainty results 


The results show that as little as 1 5% parametric uncertainty causes a divergence from stability. The 
implication is that the D1 controller is quite sensitive to parametric uncertainty. 

Next, uncertainties are applied individually to the previously mentioned four coefficients. This 
iteration, on the p-analysis procedure, shows that uncertainty in N r produces the least system stability 
sensitivity; uncertainties in Lp lead to the most sensitivity. Figure 5.15 summarizes these results. 
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Figure 5. 15 Maximum uncertainty tolerances for stability 


5.2 Linear Quadratic Robustness Analysis Applied to the X- 
38 


5. 2. 1 In troduction 

A useful and easily understood performance index, J w , is proposed by GhaouP, et al. Their 

time domain performance index is simply the value of the usual linear quadratic performance index. 
However, Ghaoui, et al., have shown that when this performance index is used for worst- case analysis 
by taking the worst initial condition vector of unit magnitude, it is a simple, yet powerful metric for 


1 . Comparing the performance of different controllers used with a given plant, 

2. Determining the worst disturbance histories for a given open- or closed- loop plant, and 

3. Determining the worst parameter changes for a given open- or closed- loop plant. 


Using this performance index, the controller synthesized in Section 4.8 is analyzed for 
performance and robustness. First, the performance of the LQG controller is compared to that of full- 
state feedback; i.e., the LQR controller. Then, the same index is used to evaluate the robustness 
properties to parametric uncertainty as well as to sensor noise and external disturbances (in particular, 
side force gust). Analysis is extended to a nonlinear system with an LQR controller. Both the control 
surface positions and rates are allowed to saturate, and the guaranteed domain of stability is obtained. 
Finally, this nonlinear analysis is further extended to a simple control surface actuator failure analysis. 


5.2.2 Performance Analysis 

Given a stable and observable linear dynamic system, 

)&=Ax+Bu, x(0)=x o , 

y =Cx 


(5.14) 
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J w is defined as 


J n 


= max 


1 > 

x'x 

A 0 A () 


This is equivalent to 


J n = max[/l(S)], 

where S is obtained by solving the steady- state Lyapunov equation 


SA + A ? S + C 7 QC = 0 


and the state initial condition, x „ , is the corresponding eigenvector 

x () =eigvec( S). 


(5.15) 


(5.16) 


(5.17) 


(5.18) 


The following index is used throughout the controller analysis: 

./„ =j^(p' +<p 2 +8a 2 + 8&+Sr 1 +S&}]t (5.19) 

The value of this index is J„ L qg = 1 -25e+5 for the designed LQG controller. If we assume lull- 
state feedback (LQR controller), the index decreases to J„lqr = 2.93e+4. Thus, the LQR controller is 
approximately four times better than the LQG controller using the performance index specified in 
Equation (5. 19). The difference in performance comes from the fact that the LQG controller has to 
estimate unmeasured states using an observer, whereas the LQR controller uses all “perfectly 
measured” states for feedback. 

The corresponding worst initial condition vector of unit length for the LQG controller is 
calculated as 


x fl = \p a A) C 0(1 &„ 8i& 8r () fig* 4" \ 

= [4.5c* — 3 -3.3c -4 - 1.7c — 4 -3.2c-4 -1.5c-6 . (5.20) 

- 5.3c -5 -1.6c -6 - 5.4c- 5 L3c — 1 9.9c-l] 

This worst initial condition vector indicates that performance is highly sensitive to the desired dynamics 
states, /4 e \ • Among the vehicle states, sideslip angle is the most sensitive state to the overall 

performance. 


5.2.3 Robustness Analysis - Parametric Uncertainties 

One application of this time-domain quadratic performance index is to measure the performance of the 
system with parameter changes. Conceptually, this process first finds the worst direction in the 
parameter space and stretches parameter variations in this direction until the system becomes unstable 
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(i.e., J n , -» oo ). A brief summary of how this robustness criterion applies to a compensated system is 
paraphrased from the original work of Ghaoui 29 , et al. Consider a linear system 

5&=Ax+Bu, x(0)=x 0 , ^ 521 

y = Cx 


with a dynamic compensator 


= A f x f +B f y 5 , x(0) = x (l , (522) 

y, =c x 

where the c subscript denotes controller values and the .v subscript denotes measured properties. Then, 
the performance index is rewritten as 


J w = max 
*<> 


JVQy^ 

X ? X 

A 0 A () 


which can be calculated easily from 


J n . - max(A(S)) 


where S is the submatrix in the solution of the steady- state Lyapunov equation 

S„A ( , + AlS a +Q„=0 

where 


A 

BC ( 

, Q = 

CQC 

0 

and S„ = 

'S 

(•)■ 

_B,C. 

A, . 


0 i 

c 7 RC f 


LO 

s,_ 


(5.23) 


(5.24) 


(5.25) 


with 


x 


a 


X 

X 

c 


and Q and R are the performance weighting matrices that define the performance index in Equation 
(5.24). The a subscript denotes augmented properties. Then, the performance of the system with 
parametric variations is obtained as follows. 

Let p be the vector of plant parameters of interest, where 

A = A(p), B = B(p), C=C(p). (5.26) 


80 



The parameter vector, p, can be broken down into combinations of the nominal value of p, Pn„m, and 
the variation from the the nominal value, Ap, as 


P = P,„„„ + ? P 

Then, a scalar measure of simultaneous changes in all parameters is defined as 

c r (p)= V( ? P) ? S" 2 ? p 


(5.27) 


(5.28) 


where 


E = diagonal matrix of standard deviation. 

Now, we can determine the Ap that maximizes J u for a specified value of o in Equation (5.28). 

Ghaoui 29 , et al., have shown that a necessary condition for the maximum is 

?P"P-P.„»=' s » (».»> 

where a is a unit vector in the direction of the gradient 


J = S 

l> 


a/J 

dp 


(5.30) 


which is evaluated at p and 


dJ K 

dp, 


is calculated as 


dJ „ 

dp, 



C t Q 


L V, 




dp, 


where P is determined by the Lyapunov equation 

AP + PA t + x„xd =0- 


(5.31) 


(5.32) 


Note that 


P = |”xx'rfr (5.33) 

This equation starts with the worst initial condition vector x» (i.e., the initial condition that maximizes ./ u ). 

Twelve parameters are selected for parametric uncertainty analysis. They form the vector of 
plant parameters, p, in Equation (5.27) as 

P = h y s L, i„ L, L„ L, A', N„ N, N„ nJ (5.34) 


81 



Then, the value of three standard deviations for each parametric uncertainty is assumed to be equal to 
the magnitude of the nominal value of each parameter. 

For the designed LQG controller, the gradient of the performance index with respect to these 
1 2 parameters was calculated as 

— S=[-2.8el 3.3 -2.2e7 1.8e5 6.6e2 -5.7e4 

dp L • (5.35) 

1.2e4 - 1.5e3 3.0e2 5.4e3 -3.5e4 2.9e4] 

The magnitude of this gradient divided by the nominal ./„ is 1 75, indicating that the performance index 
increases 175 times for a one-sigma change in the worst direction for the parameter space. 

Figure 5. 16 shows \/J» versus a, where a is defined in Equation (5.28) for the system with the 
LQG controller. 



sigma 


Figure 5.16 1/J W versus a for worst parameter change 


System instability occurs at a = 1 .65 in the worst direction in the parameter space. At the stability 
parameter margin, i.e., a = 1 .65, 

Ap r S' 1 = [- 2. 1 e - 6 2.5e - 7 -1.7 1.3e-2 5.0e-5 -4.3e-3 ( 

9.0e -4 -l.le-4 2.3e-5 4.1e-4 -2.7e-3 2.2e-3] 


82 




This indicates that the most important major contribution is from a decrease in Lp and, to a 
lesser degree, a decrease in L,,. The sensitivities of performance due to all other parameters are 
negligible as compared to these two parameters. 

5.2.4 Robustness Analysis - Disturbance 

So far, all perfonnance analyses assume no disturbances are present. However, the presence 
of noise is inevitable in physical systems. In this section, the original performance criterion is modified so 
that both process noise, vv. and measurement noise, v, are addressed in the performance index. To 
achieve this, the worst disturbances are assumed to be feedbacks of the augmented state, x a , where the 
gain matrices are determined by the solution of the Riccati equation 


s„a„+a:s„+q„+-s„grXs, ; =o, 

n 


(5.37) 


where S a , A a , and Q a are defined in Equation (5.25) and 
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G 0 

, R = 
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0 “ 

0 B 

’ a 
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R, 


M’ = K . 


v = K, 


with 


K, 


X„ 

K,; 

1 

R 'G 7 0 

K 


x„ 

K w . 


o r;'b , 7 


Now, J w is given by 


,/ u = max [A(s)]+/ ? *r, 


where I'l is a Lagrange multiplier and W is determined from 

ir = //(l'(B;Nk , 1]" 

and P„ is the solution to the Lyapunov equation 

(A „ + G„ K„ )P„ + P, , (A, + G , K a ) r + x„ 0 x:, 0 = 0 


(5.38) 


(5.39) 


(5.40) 


(5.41) 


(5.42) 


and 


K 


a 



(5.43) 


First, the original system defined in Equation (4.3 1 ) is modified to accommodate a disturbance 
due to side force gust. 
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( 5 . 44 ) 


As 


^ disfidist 


B rf ,. s ,U rf) , ; + G d,s,™dis< 


where 


^ dm ® an}; ’ ^ dist ^auy 1 ^ f/jw ®o«ir 

G (/m - [- ^ /V T -Lp -N p 000000 O] 7 (5.45) 

W rf/Sf = P 

The difference from the original system is the addition of the v term. This formulation is 

interpreted as side force translated into the sideslip angle, and this disturbed sideslip angle acts as a 
control although it is trying to destabilize the system. In our example, we select max I (3 I = 1 degree for 
transonic flight. 

Next, another modification to the original system is made to accomodate the measurement 
noise. Adding sensor errors to continuously changing scale factors and biases to the measurements 
vector is the technique used. 
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(5.46) 


Also, we assume the worst disturbances are feedbacks of the augmented state x a , where gain matrices 
are determined by the solution of the Riccati equation specified in Equation (5.37). In our example, R w 
= 1 and R v = diag(3,l ) are selected and the corresponding value of rj is 6.84 x 1 0 4 by interpolation. 
Disturbances are given by the positive feedbacks defined in Equation (5.38); and the feedback sideslip 
disturbance vector, K„ , and sensor noise matrix, K,, are calculated as follows: 

K, = [- 5.9c -3 2.9c - 3 2.1c-2 -1.0c-4 6.8c-6 2.5c - 4 ~1.4c-5 -5.1c-4 -l.8e-4 -!.8c-2 

-2.V- 3 - 1 ,2c - 3 - 1.2c- 2 1.3c -3 -1.3c-6 -4.2c -5 9.9c-6 3.6c-4 7.8c-4 2.0c-2] 

(5.47) 

r2.0c-5 - 6.7c- 6 -1.8c- 5 -l.9c-6 -2.5c-8 -9.1c-7 -9.1c-9 -3.2c- 7 -6.3c-6 -1.8c-5 

K = 

' [].4c-5 4.9c-5 1.5c- 3 -8.lc-5 -2.7c- 7 -1.0c- 5 -l.8c-6 -6.8c-5 8.8c-5 -8.6c-4 

- 5.0c- 6 - 1.1c- 7 -5.4c -6 -2.5e-6 -1.3c-9 -6.8c-8 8.6c-9 3.0c-7 3.3c-6 3.2c -6‘ 

- 3.8<' - 4 -8.6c- 6 -8.8c- 4 6.5c-5 2.8c-7 1.0c-5 l.lc-6 4.lc-5 ~5.5c-5 1 ,2c- 3 


The value of increases to 1.69e+5, a 35% increase from the nominal case value of 1.25e+5. The 
worst unit initial condition vector is 

x = [4.8e — 3 - 2 ,9e — 4 -2.6e-4 -3.4e-4 -1.3e-6 

17 • (5.48) 

— 4.7e—5 -\3e-6 -4.4c -5 1.3e-l 9.9e-l 


The magnitude of the initial measurement noise is obtained by substituting Equations (5.47) and (5.48) 
into Equation (5.38). The ratio between the actual values to the noise for sideslip angle and bank angle 
is calculated as 5.4%, and 17.4%, respectively, at t = 0. 
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5.2.5 Domain of Stability for the System with Actuator Saturation 

Although the time domain criterion we have chosen is a useful tool for controller robustness 
analysis, its major drawback is that the entire structure is defined using a linear system assumption. 
However, most aerospace systems are not linear; so this point is especially significant since the 
performance of a DI controller is sensitiv e to the available control power. This concern essentially 
arises because the control inputs, u(x), is proportional to the inverse of the magnitude of the control 
distribution function. 

Therefore, the control surface position and rate should be included when the boundary of the 
stability region is considered. Although this is a significant weakness of DI methodology, relatively little 
research has come to our attention 10 " 3 '. Recently, Tarbourieclf”, et al., published a technique to 
compute a guaranteed domain of stability for a system subject to position- and rate- limited system 
inputs. This study shows two different approaches: the algebraic Ricatti equation ( ARE) approach and 
the linear matrix inequality (LMI) approach. Here, the simpler and more widely used approach of the 
two, the ARE approach, is used to find the domain of the system stability. 

The objective is to find the largest possible Lyapunov stability parameter, p , such that the 
closed-loop system is locally stable in the largest Lyapunov level set. In other words, we seek the 
largest domain of initial state vectors that produces a stable solution to the given Ricatti equation defined 
by Tarbouriech 11 , et al. Though this methodology produces an optimal solution by solving a given 
Ricatti equation, this optimal solution depends on the choice of state and control weight matrices, which 
is always the case for Ricatti solutions. Therefore, the largest domain of initial condition predicted by 
this approach does not grantee the largest stable initial condition domain globally, and the computed 
stability domain may still be a conservative prediction. 

Actuator position and rate limits were previously listed in Section 4. 1 .2. Using the DI 
formulation, the actuator dynamics are not accessible directly. Instead, the desired yaw and roll 
acceleration and their time derivatives are used to limit the control inputs. For the yaw axis, the desired 
yaw acceleration is bounded by 


\tSt 


N 


da 


+ 


A^-I^L = 0.485(/W / sec 2 ) 


(5.49) 


and the time derivative of the yaw acceleration is limited to 


gt 


max 


+ A aM,,x = 0.990 (rad! sec 3 ). 


Similarly, for the roll axis. 


m L = L &Nmox +L A 5r Ls = 2 -69 (rod/ sec 2 ) 


(5.50) 


(5.51) 


and 


fit" 


= L 


6a 


+ L 


,|<54 max = 5.66 (rad! sec 1 ) 


(5.52) 


Assuming full-state feedback with the regulator gain in Equation (4.39) and no parametric uncertainties 
and disturbances, the maximum value of the stability parameter (p nm =) is 14.5. When we consider the 
domain of stability in the two different states - bank and sideslip angles - simultaneously, a 3-D plot is 
obtained (Figure 5.17). The system is stable up toa 48.1° bank angle assuming no sideslip angle. 
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Similarly, the system is stable up to 1 .59° sideslip assuming no bank angle. The system is not 
guaranteed to be stable outside of this domain. 



Sideslip Angle (deg) u u Bank Angle (deg) 


Figure 5.17 Stability boundary 


5.2.6 Change in Domain of Stability due to Control Surface 
Actuator Failure 

This nonlinear analysis next examines the stability domain in the event of a control surface 
actuator failure. This process changes the limits on the position and rates of the input vectors. Six 
failure modes are considerd, and results are compared to the nominal case where no control surface 
actuator failure has taken place. The nominal and failure mode cases that are considered are the: ( 1 ) 
nominal case, (2) one aileron failure case, (3) one rudder failure case, (4) one aileron and one rudder 
failure case, (5) two aileron actuators failure case, and (6) two rudder actuators failure case. Also, as in 
the previous section, perturbations on the inital conditions are limited to bank and sideslip angles. Figure 
5.18 summarizes the results and their associated cross-sectional top view. This figure clearly illustrates 
the difference in stability Lyapunov levels due to different constraints on control inputs. As predicted, 
the nominal case, where none of the control surface actuators are failed (and the largest amount of 
control power is available), shows the highest level of stability. The ratio, compared to the nominal 
stability domain, for all six cases is summarized in Figure 5.19. 
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Figure 5.19 Area of stability comparison due to actuator failure 
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6 Theoretical Foundations 


6.1 Basic Forms of Dynamic Inversion 

Dang Vu asserts that the essentials of the DI approach are most easily understood in terms of an 
SISO system. Since we followed that approach in our learning process, as described in Section 3, we 
have postulated a class of linear systems affine 1 in control that is represented by the following 
mathematical form: 


&=f(x) + g(x)u (6.1) 

y = h(x) (6.2) 

where /(.v) and g(.v ) are smooth vector fields on R" and h(x ) is a smooth mapping function that maps R" 
into R. This system is feedback linearizable of relative degree r if state and input transformations exist. 

z = <D(x) zeR r 

y = a( x) + P(x)\/ i/gR 

where (3(x) * Oand 0(.v) is a diffeomorphism" that transforms Equation (6.1) into a controllable linear 
system. 


&=Az + Bv 

Following Dang Vu’s rather succinct development (with some clarifications), we differentiate the 
nonlinear output equation (Equation (6.2)) with respect to time and obtain 

If the coefficient of it is zero, we continue with successive differentiations in the same fashion until a 
nonzero process coefficient surfaces. Dang Vu’s succinct notation, which uses the Lie derivative from 
differential geometry for these repeated differentiations, is useful. 


'“Affine in control” means that all transformations of finite system parameters remain finite under the 
controller’s action. 

"A diffeomorphism is said to occur when the scalar components of a mapping (or transformation) F are r 
times differentiable with r > 1 with respect to the scalar components of x (with a discrete time mapping Xk 
= F Xk+i ). The mapping must also be invertible; that is, x k+1 = F ' 1 x k+1 must hold. Invertibility implies that 
F -1 exists. The scalar components of the inverse must likewise be r times differentiable. 34 
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(6.3) 


Equation (6.3) is the Lie derivative of the scalar function h with respect to the vector field /. Higher- 
order derivatives have a similar fomi. 

L k f h(x) = L f ^L k f ~ 1 h(x)^ (6.4) 

Using this Lie derivative notation, the output equation can be rewritten as 

^^L&=^f(x) + g(x)u] = Lfh(x) + L g h(x)u. (6.5) 

OX oX L 

If the second derivative in Equation (6.5) is 0 - that is, Lgh(x) = 0 - a second differentiation yields 

flt=Ufh(x) + LgLfh(x)u. (6.6) 

The differentiations end when Lgdf 1 h(x) = 0 fork=l, ... , r- 1, but Lgllf 1 h(x)^0. The last 
derivative of the sequence for the output response is then 

y^ r> = Lfh(x) + LgLf~ 1 h(x)u . (6.7) 

The number r is called the relative degree of the original control equation (Equation (6. 1 )). 

If we define our coordinate transformation in terms of the Lie derivative, 

z k = Q> k (x)=L k ~ 1 h(x) for A' = /, 2, ..., /•, 

the resulting transformed set of equations is linear, of dimension r, and in a companion form called the 
Brunovsky canonical form. 

'0 10 L 0] \0' 

0 0 1 L 0 0 

&= 0 0 L 1 0z+0v = Az + Bv (6.8) 

0 0 L 0 1 M 

0 0 L 0 0] [l_ 

where v = Llf h(x) + LgL r f 1 h(x)u . Vu points out that exact linearization is possible when the relative 

degree r is equal to the order of the system n and the lineanzed system (Equation (6.8)) is both 
controllable and observable. 

Since we are interested in obtaining for the linearized system a control law that will impose 
desired behavior on the original nonlinear system, it is also necessary that we carefully examine any 
control law in terms of stability of the nonlinear system. Obviously, that requirement suggests that we 
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transform any postulated control law back into the original coordinates. Symbolically, we can write this 
transformation as 


u = a( x) +(3(x) v 


with 


a(x) = - 


Lgl^f 1 h(x) 


and P(x ) = 


1 _ 

LgUfMx) 


6.2 Stability and Robustness Analyses 

In this section, stability and robustness analyses are described for the DI controller. The most 
commonly used methodology used to analyze robustness of linear systems using DI controllers is based 
on the structured singular value (|i) and a technique now widely described in the literature as (i-analysis. 
However, p-analysis is not the only method used to examine the stability and robustness of Dl-based 
controllers. 


6. 2. 1 Linear Fractional Transformations 

LFTs are used to integrate parameter variations (uncertainty) into the system under p-analysis. 
As the first step in the p-analysis procedure, all parameter variations are collected into an uncertainty 
matrix in LFT form. Then, p- analysis looks for the smallest variation in these parameters that drive the 
system to instability. In this section, the methodology to create an LFT form is explained. 

Suppose we have a linear system that is described by the following: 


#= ax + bu 
y = x 


( 6 . 9 ) 


Now, we assume that the value of a varies between a and a + . 


a < a < a + 


( 6 . 10 ) 


where 


a' = lowerlimit of variation in a 
a + = upper limit of variation in a 

This relation can be rewritten in terms of nominal value of a , a nom , as 


a = a nom + 


kA 


1-k 2 h 


( 6 . 11 ) 


where 



a nom = nominal value of a 


ki = 


2(a + -aj(a-a j 


a + - a" 


ko = 


|a + -t-a j -2a 


a + - a' 


- 1 < 5 a <1 


When the upper and lower variations from the nominal values are equal - that is, when 


a + - a nom 


a - a 


nom 


- then the previous equation simplifies to 


a = a 


nom 


+ /c8 s 


(6.12) 


where 


a nom = nominal value of a 
k = a + - a nom =a nom -a~ 

~1<h a <1 

The perturbation in a described in Equation (6.1 1 ) is integrated into Equation (6.9) and is now 
expressed in state- space form by introducing the fictitious terms Z a and w a as 




a 

k i 

b~ 

X ” 


- 

1 

k 2 

0 

w a 

J _ 


1 

0 

0 

u 


The block diagram in Figure 6. 1 (above) corresponds to the state-space expression (Equation 
(6. 13)). Here, we have completely separated what is known, the G matrix, from what is uncertain, A. 
We treat A as uncertain; but. we do know its range: -1 < A < 1 . For now, uncertainty in a alone is 
considered. The next case to be considered is the case when there is an uncertainty in h . This case, 
variation in b , is expressed similar to the previous case as 




a 

k 'i 

b 


Zb 

- 

0 

*2 

1 

w b 

y _ 


1 

0 

0 

j u 


where 


(6.14) 
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b + -/b)(<b-/r 

b + -b~ 


ib + + b~)-2b 





Figure 6. 1 Linear Fractional Transformation block diagram 


Finally, combining the expressions for uncertainty in a and in b , the following state- space 
expression results. The accompanying block diagram (Figure 6.2) is a companion to Figure 6. 1 . 


'# 


a 

k l 

k l 

b 

X 



1 

k 2 

0 

0 

w a 

Zb 


0 

0 

k 2 

0 

w b 

_y . 


1 

0 

0 

0 

u 


(6.15) 



Figure 6.2 Companion to the Linear Fractional 
Transformation block diagram 


Now, A is no longer a scalar uncertainty but is a 2 x2 diagonal matrix with normalized 
uncertainty terms. 
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(6.16) 



0 

§6 


Again, we have separated what is known from what is uncertain but bounded. Since A is no longer a 
scalar, we have to make another choice for a norm: the maximum singular value. It is not difficult to see 
that (7 (A) < 1. A very important observation is that the uncertain element A has a fixed structure: a 
diagonal matrix consisting of the individual uncertainties in a and b. Thus, unstructured uncertainty at 
the component level has become structured uncertainty at the system level. LFTs are the mathematical 
tools that allow us to provide this systemic structure for the uncertainty. 

6.2.2 Other Types of Uncertainty Models 

Other than the parametric uncertainty described in Section 6.2. 1 , there are at least two other 
types of uncertainty models. These are Unmodeled Dynamics or Uncertainty at the Input and 
Uncertainty at the Output models. The major difference between parametric uncertainty and the other 
two types of uncertainties is that parametric uncertainties are real- valued while the other types of 
uncertainties are complex- valued perturbations. The two types of uncertainty models are explained in 
the following subsections. 


6. 2.2. 1 Unmodeled Dynamics (Uncertainty at the Input) 

During the linearization process, higher-order terms in aircraft equations of motion are ignored. 
Also, other uncertainties arise due to aeroelasticity, control surface variations, and vehicle flexibility. 
Usually, the plant model is a good system representation tenn at low- to mid- frequency inputs, but 
modeling uncertainties become larger with high-frequency inputs. Instead of attempting to include all 
modeling uncertainties, the modeling uncertainties are treated as additives to the plant inputs. Figure 6.3 
shows a general block diagram representing this approach to accounting for uncertainty at the plant 
input. 



Figure 6.3 Unmodeled Dynamics block diagram 


This type of uncertainty is parameterized with two elements, W m and A ,,,. W m is a weighting 
transfer function (assumed to be known) that reflects the amount of uncertainty in a model with respect 
to frequency. The other parameter, A„„ is a stable unknown transfer function that nevertheless satisfies 
the condition 1 1 A J I „ < 1 . 


6. 2. 2. 2 Uncertainty at the Output 

Similar to Unmodeled Uncertainty, uncertainty of the measurements is modeled as Uncertainty 
at the Output. A block diagram representing how this type of uncertainty is modeled is shown in Figure 


6.4. 






Figure 6.4 Uncertainty at the Output block diagram 


6.2.3 Structured Singular-Value Analysis (p- Analysis ) 

Now that we know how to represent uncertainties in the system using LFTs, we must turn our 
attention to analyzing the robustness of systems modeled in this fashion. We follow the most common 
practice today by basing our analysis on the structured singular value, |i, and using available software 
tools to perform the p- analysis. The technique is based on the following theorem: 

Robust Stability <=> p A ( jt o)) <1 V to (6. 1 7) 


where M n is the left upper comer block of A7; i.e.. 


M{j co) 


Mu{j(o) M 12 {ju>) 
M 2 i{j(o) M 22 (joi) 


(6.18) 


and the function p A is defined as 

( M ) = min{o( A): Ae A,ctef (/-MA) = 0} (6 ‘ 19) 

where A = j diag ( A-/ , A 2 , A n )]. 

According to this theorem, p A is a function of M that depends on the structure of A . p A is the 
reciprocal of the smallest A (where we use o as the norm) we can find for the set A that makes the 
matrix / - MA singular. If no such A exists, |i A is taken to be zero. 

The general framework and the |i- analysis transformation are shown in Figure 6.5. 
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Figure 6.5 General framework and \L-analysis transformation 


Even though the function ji A is defined, we must still calculate it. Unfortunately, no exact 

calculation algorithms exist. So, we must calculate its upper and lower bounds. Normally, the upper 
bound is used since these values of |i are “safer” (that is, they are more conservative). The upper 

bound is defined as 

Ma (M)< inf g{dMD^). (6.20) 

- DeD ' ' 

where D is the scaling matrix. Figure 6.6 geometrically illustrates the effect of D -scales. 



Figure 6. 6 The effect of D -scales 


Another important feature of the upper bound is that it can be combined with the FL controller 
synthesis technique to yield a u- synthesis method. Note that the upper bound, when applied to transfer 
functions and maximized across frequencies, is simply a scaled H. norm. 

The steps needed to test the robust stability using p -analysis are as follows: 

1 . Construct the interconnection structure, M, which is a known linear system. 

2. Define a structured perturbation set, A. 

3. Combine M and A to form the feedback system shown in Figure 6.5. 

4. Calculate a frequency response of M. 

5. Calculate the upper and lower bounds for p. 

6. Find the upper bound peak value. 
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7. If (V«* < 1 : pass; if \i, Ka k > 1 : fail. 
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